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Preface 



This volume is a collection of original and review articles which reflects many 
developments on the cutting edge of research in operator theory and its appli- 
cations, and related areas of mathematics: multivariable operator theory and its 
connections with multi-dimensional systems and the function theory of several 
complex variables, operator theory in indefinite inner product spaces, distributed 
parameter systems, mathematical physics (inverse scattering, Dirac-type equa- 
tions, Schrodinger operators), wavelets (including connections with the representa- 
tion theory for operator algebras), and numerical analysis. Review articles include 
a report on recent achievements and future directions of research in the area of 
operator theory and its diverse applications. 

The collection appears as proceedings of the International Workshop on Op- 
erator Theory and its Applications (IWOTA) which was held at Virginia Tech, 
August 2002, and reflects lectures, talks, and discussions which took place at the 
workshop. 

We gratefully acknowledge the following agencies for their financial support 
of IWOTA 2002: 

NSF (National Science Foundation) 

Department of Mathematics, Virginia Tech 
College of Arts and Sciences, Virginia Tech 
Research Division, Virginia Tech 
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IWOTA 2002 and 

Recent Achievements and New Directions 
in Operator Theory and Applications 

Joseph A. Ball, J. William Helton, Martin Klaus, and Leiba Rodman 



Abstract. This is a report on IWOTA 2002 and recent achievements, chal- 
lenges and directions for future research in operator theory and its applica- 
tions, as reflected at the International Workshop on Operator Theory and 
its Applications (IWOTA), held at Virginia Tech, Blacksburg, Virginia, from 
August 6 till August 9, 2002. The report incorporates input received from 
workshop participants and others whose fields of research are included in the 
IWOTA 2002 themes. We thank all who helped us to prepare this report. 



1. The workshop IWOTA 2002 

The primary objective of IWOTAs, which have met biannually since 1981 and 
more frequently in recent years, is to bring together major researchers in the area 
of operator theory and related fields, including applications in engineering and 
mathematical physics. IWOTAs emphasize cross-disciplinary interaction between 
mathematicians, engineers, physicists, and scientists in other fields in which results 
and methods of operator theory have significant applications. IWOTA 2002 had 
14 plenary lectures; the plenary speakers and the titles of their presentations were 
as follows: 

Damir Arov, The class of strongly regular J -inner matrix functions 
Vladimir Bolotnikov, On interpolation problems for contractive-valued func- 
tions on certain domains in C n 

Stephen Boyd, Optimization over linear matrix inequalities 
Sarah Ferguson, A new look at H 2 of the bidisk 
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torization, matrix functions, linear matrix inequalities, distributed parameter systems, direct 
scattering, inverse scattering. 
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Tryphon Georgiou, High resolution spectral analysis: advances and applica- 
tions 

Israel Gohberg, Orthogonal systems and convolution equations 
Palle E.T. Jorgensen, Factorizations and decompositions for operators moti- 
vated by wavelet analysis 

Marinus A. Kaashoek, Two new additions to the commutant lifting theorem 
Vadim Olshevsky, Efficient algorithms for structured matrices 
James Rovnyak, Spectral problems for some indefinite generalizations of can- 
onical differential and difference equations 
Lev Sakhnovich, Matrix finite- zone Dirac type equations 
Victor Vinnikov, Systems on Lie groups and nonself adjoint representations 
of Lie algebras 

Alexander Volberg, Electric intensity capacity and operator theory 
Ricardo Weder, The time- dependent approach to inverse scattering 

In addition there were 73 talks grouped by special sessions: Interpolation 
theory (V. Bolotnikov, organizer); Scattering and Inverse Scattering (M. Klaus, 
organizer); Operator Theory and Applications (H. Langer, organizer); Wavelets 
(D. Larson, organizer); Structured Matrices (V. Olshevsky, organizer); Distributed 
Parameter Systems/Operator Theory (O. Staffans, organizer); Multivariable Oper- 
ator Theory (V. Vinnikov, organizer). The speakers and titles of their presentations 
were as follows: 

Interpolation Theory: 

Yury Arlinskii, On von Neumann’s problem in the extension theory of non- 
negative operators 

Animikh Biswas, Intertwining lifting of operators 

Tiberiu Constantinescu, Szego kernels and orthogonal polynomials in several 
noncommuting variables 

Aad Dijksma, Nudel’man’s problem in an indefinite setting 
Harry Dym, Application of Riccati equations to singular interpolation prob- 
lems 

Gilbert J. Groenewald, J -spectral factorization for arbitrary rational matrix 
functions 

H. Turgay Kaptanoglu, Radial differential operators and Bergman projections 
on diagonal Besov spaces on the ball 

Alexander Kheifets, Defect and equality in boundary interpolation 

Linda Patton, Pick interpolation on the boundary of the bidisc 

Dirk Pik, The Kalman- Yakubovich- Popov inequality and infinite- dimensional 

discrete-time systems 

Andre C.M. Ran, A fixed point theorem in partially ordered sets and some 
applications to matrix equations 

Eduard Tsekanovskii, Realization theory of Herglotz-Nevanlinna matrix- 
valued functions, interpolation problems and explicit system solutions 
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xi 



Scattering and Inverse Scattering: 

Tuncay Aktosun, Inverse scattering with partial data 

Jared Bronski, Eigenvalue asymptotics , passive scalar intermittency and small 
ball estimates for fractional Brownian motions 

Ira Herbst, Quantum scattering for a class of long-range magnetic fields at 
high energies 

Peter Miller, Scattering theory for nonself adjoint Zakharov- Shabat operators 
in the semiclassical limit: some facts arid open questions 
Alexander Rybkin, Higher order WKB-type phase for the Schrodinger equa- 
tion in ID and scattering on long-range potentials 
Ken Shaw, Eigenvalues and pulse shapes for Zakharov- Shabat systems 
Alexey Tikhonov, Functional model for operators with spectrum on a curve 
Vadim Tkachenko, 1-d periodic differential operator of order 4 
Alexander Tovis, Semiclassical limit of a focusing nonlinear Schrodinger equa- 
tion and related scattering problems 

Cor van der Mee, State space methods for solving inverse scattering problems 
on the line 

Rudi Weikard, On the inverse resonance problem 

Peter Yuditskii, Functional models for almost periodic Jacobi matrices 

Operator Theory and Applications: 

Miron Bekker, Helson-Szego theorem for operator-valued weights 

Aad Dijksma, The Schur algorithm for generalized Schur functions and J- 

unitary matrix polynomials 

Seppo Hassi, A reproducing kernel space model for N K - functions 
Lazio Kerchy, Generalized Toeplitz operators 

Heinz Langer, Self-adjoint block operator matrices with non-separated diago- 
nal entries and their Schur complements 

Matthias Langer, Variational principles for eigenvalues of pencils of un- 
bounded operators 

Victor Lomonosov, The Bishop-Phelps theorem fails for uniform non-self- 
adjoint dual operator algebras 

Annemarie Luger, A characterization of generalized poles 

Konstantin A. Makarov, Variation of spectral subspaces under off-diagonal 

perturbations 

Dmitry J. Nicolsky, Approximation of the spectrum of a quasi-periodic func- 
tion sampled at a finite number of points distributed uniformly 
Artem D. Pulemyotov, On the generalized joint eigenvector expansion for 
commuting normal operators 

Peter Stollman, Quasicrystals and aperiodic order: ergodic properties and 
operator theoretic consequences 

Christ iane Tretter, On some quadratic operator pencils associated with 
damped dynamical systems 
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Wavelets: 

Radu Balan, The C* -algebra of nonexpansive operators and application to a 
frame redundancy result 

Yuanan Diao, The s- elementary frame wavelets are path connected 

Eugen J. Ionascu, Wavelet sets and their induced isomorphisms on [0, 1) 

David Larson, Wavelets, frames and operator theory 

Kerri Kornelson, Ellipsoidal frames 

Steen Pedersen, Exponential bases for Cantor measures 

David Walnut, Sampling and the local Pompeiu problem 

Eric Weber, Geometry of affine frames 

Structured Matrices: 

Albrecht Bottcher, The product of a structured matrix and a random vector 
Tiberiu Const ant inescu, Displacement structure, quantum channels, and ten- 
sor algebras 

Mark Embree, Pseudospectra of piecewise continuous Toeplitz matrices 
Dario Fasino, Fast orthogonalization of Cauchy-like matrices and rational 
functions 

Georg Heinig, Fast direct algorithms for Toeplitz least squares problems 
Christian Mehl, Normal matrices in indefinite inner products 
Karla Rost, Toeplitz- plus- Hankel matrices 

Michael Stewart, Orthogonalization of structured matrices using inner prod- 
ucts 

Distributed Parameter Systems: 

Birgit Jacob, When is an observation operator admissible? 

A. Karelin, On a reduction of singular integral operators with a linear-fract- 
ional involution to matrix characteristic operators 

Gabriel Prajitura, Classes of operators with fixed local spectral properties 
David Russell, Some wave and buckling phenomena in nonlinear elastic beams 
with supporting constraints 

Olof J. Staffans, Conservative continuous time systems out of thin air 
Ping Wang, The application of the Moore-Penrose generalized inverse to an 
operator equation 

Multivariable Operator Theory: 

Raul Curto, Algebraic varieties arising in truncated complex moment prob- 
lems 

Ronald G. Douglas, Quotient modules and nilpotent actions 
Michael Dritschel, A completely positive approach to Ando’s theorem 
Chen Dubi, Realization of rational functions in C N 

Xiang Fang, Dimension, multiplicity, and some other invariants for operator 
tuples 

Lawrence Fialkow, Truncated planar moment problems: the algebraic variety 
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D. Kafyuzhniy-Verbovetskii, On the Bessmertny! class of homogeneous posi- 
tive analytic functions of several variables 

Mircea Martin, Uniform approximation by solutions of elliptic equations and 

seminormality in higher dimensions 

John E. McCarthy, Operator theory on varieties 

Scott McCullough, Noncommutative polynomials and Linear- Matrix- Inequal- 
ity representations 

Gelu Popescu, Multivariable Nehari problem and interpolation 
Mihai Putinar, Privilege on strictly convex domains of C n 
Jim Solazzo, Reflexivity, factorization, and Hankel operators 
Tavan T. Trent, Corona theorem for Dirichlet space 



2. Interpolation and commutant lifting theory 

Recent achievements. Interpolation theory played a defining role in the develop- 
ment of i/°°-control theory beginning in the 1980s. As part of these developments, 
we include the Adamjan-Arov-Krein extension of the Nehari theorem (character- 
izing infinity-norm approximation of a given L°° function on the unit circle by 
a meromorphic function with prescribed number of poles inside the disk) and its 
later implementation in state-space form, as well as the development of an array 
of approaches (via commutant lifting, projective geometry/ indefinite Beurling-Lax 
theorem, reproducing kernel Hilbert spaces, unitary extensions of isometry, “fun- 
damental matrix identities” related to positive kernels, state space implementation 
via directional zero/pole and/or factorization analysis) to the solution of matri- 
cial Nevanlinna-Pick interpolation problems. A more recent advance has been the 
parametrization of the set of all solutions of a Nevanlinna-Pick interpolation prob- 
lem also meeting a low McMillan-degree constraint, originally through application 
of techniques from nonlinear functional analysis (degree theory); this result has 
important implications in speech synthesis applications. More recently there have 
been extensions of the matricial Nevanlinna-Pick interpolation theory to a context 
of time-varying systems (both in discrete and continuous time) and there have 
been various refinements of the commutant lifting theorem (more flexibility in the 
formulation of the norm constraint and in the intertwining hypothesis) which have 
led to new types of applications. In another direction, interpolation theory and 
related problems have been studied in classes of functions whose Caratheodory- 
Pick matrices have a bounded number of negative eigenvalues (with no a priori 
regularity hypotheses). 

Future directions. A newer emerging area for applications of interpolation prob- 
lems is signal processing; specifically, problems arising in the construction of an- 
tennae arrays call for the solution of some nonstandard interpolation problems. A 
general theme for future research is the characterization of interpolants to certain 
non-standard interpolation problems. In this context, as well as in the context 
of more traditional analytic interpolation theory, it is important to quantify the 
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size of families of interpolants as it relates to modeling uncertainty and resolu- 
tion in engineering applications, as well as to study the effect of interpolation 
nodes/constraints (which relate to, e.g., antennae array geometry and sensor dy- 
namics) to the size of such families of interpolants. 



3. Multivariable operator theory, several complex variables and 
multi-dimensional system theory 

Recent achievements. There has been a recent flurry of activity extending the 
interpolation and commutant lifting theory described in the previous section to 
multivariable settings. We now have a theory of Nevanlinna-Pick interpolation 
(including the case of interpolation at points on the boundary of the domain) 
analogous to the classical theory in several multivariable settings, including the 
polydisk, the ball and the general setting of a reproducing kernel Hilbert space 
with “Pick kernel” . Commutant lifting and Toeplitz corona theorems (including 
some progress on the solution of the Carleson corona theorem on domains beyond 
the disk case) have now been obtained for both the ball and polydisk setting. 
More broadly, this work includes connections with multi-dimensional system the- 
ory, model theory for tuples of commuting operators in various classes and more 
general multivariable settings for Lax-Phillips scattering. The solution of the inter- 
polation problem (or the characteristic operator function or the scattering opera- 
tor, depending on the context) is usually realized explicitly as the transfer function 
of some sort of conservative, multi-dimensional linear system. There now has ap- 
peared some work which extends these ideas beyond a specific domain (such as the 
polydisk or the ball) to a general domain in C n defined as the points of contrac- 
tivity for a given matrix multivariable polynomial; the classical Cartan domains 
of the first three types arise as special cases. Related work deals with a noncom- 
mutative version of Schur-class of functions (the unit ball of the “Cuntz- Toeplitz 
algebra”), the elements of which are represented as formal power series in some 
number of noncommuting indeterminants, multiplication by one of which gives rise 
to a contraction operator on the Fock space commuting with the right creation 
operators. Results in the noncommutative theory can be symmetrized to produce 
the corresponding results for the multiplier class of functions on the ball studied 
by Arveson. More abstract approaches are also being pursued for the commutative 
case which use the language of modules over function rings and homological alge- 
bra. There have also been recent breakthroughs in the understanding of the various 
types of multivariable classes of functions of bounded mean oscillation (BMO) and 
in estimating the norms of big and little multivariable Hankel operators, as well 
as continuing progress on various sorts of multivariable moment problems (which 
have important applications in signal processing). 

Future directions. There already exist a number of older results (due to Cartan 
and others) on holomorphic extension of a given complex- valued function defined 
on a subvariety to a larger domain which were obtained by homological-algebra 
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methods. An important direction for future research is the assimilation of these 
older results based on homology methodology into the operator/system theory 
framework which has led to the recent analogues of Nevanlinna-Pick interpolation 
discussed above. 

There already exists a now relatively old result due to Ando on the existence 
of a commuting unitary dilation for a commuting pair of dilations. We now appear 
to be on the verge (from a number of independent directions) of getting an explicit 
understanding of the geometric structure of this dilation. 

We expect areas of continuing activity to be the development of classical 
interpolation theory on nonclassical domains. More broadly, we expect to see the 
integration of operator model theory and system theory in more general contexts. 
Currently, we have such an integration for polydisks and interpolation in the Schur- 
Agler class (a class somewhat smaller than the Schur-class of functions mapping 
the polydisk in C d ( d > 2) to the unit disk in C), for balls and the associated 
Arveson class of contractive multipliers with respect to the kernel kd{z,w) = 1/(1 — 
(z,w)), and for the class of contractive bundle maps between two Hardy spaces 
consisting of holomorphic sections of parahermitian vector bundles over a real 
Riemann surface. There have been generalized settings of operator model theory 
due to Curto, Vasilescu and others for which the connections with system theory 
and Lax-Phillips scattering theory are not clear. 

By taking the Riemann surface to be the double of a finitely-connected planar 
domain, this last setting includes the case of Hardy spaces over planar domains 
studied by Abrahamse and Douglas in the 1970s. Since the original work of Abra- 
hamse on the Nevanlinna-Pick interpolation problem on a finitely-connected planar 
domain, there has been some but limited progress in understanding the structure of 
the set of all solutions, or in obtaining the result via system-theory ideas (whereby 
one arrives at the solution realized as the transfer function of some linear system). 
A general setting for these issues is the Nevanlinna-Pick interpolation problem for 
contractive multipliers on a reproducing kernel Hilbert space, where the kernel is 
not a “Pick kernel” for which the Nevanlinna-Pick interpolation theorem holds 
in essentially the classical form. We expect future work eventually to unveil the 
mysteries currently surrounding this area. 

While there has been much work in the multivariable-operator-theory liter- 
ature on commuting contractions or commuting nonselfadjoint operators close to 
selfadjoint, this work has not made any connections with the extensive literature 
on hyponormal operators, which until now has used different techniques and dealt 
with different issues. The theory of hyponormal operators somehow lies deeper - 
more work needs to be done here to make connections with the recent developments 
in other parts of multivariable operator theory. 

We now have a couple of isolated instances of multivariable commutant lift- 
ing - this should be understood in a broader, more fundamental setting; in this 
connection, the Taylor functional calculus for operator tuples needs to be bet- 
ter known and more broadly used. This general commutant lifting theorem ap- 
pears now to be a key piece (along with the 9-technique) toward an operator- 
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theoretic simplification and proof of the Carleson corona theorem on more general 
domains. 

We also expect a more complete theory of multivariable moment problems 
to emerge. 

For all these directions we expect to see a continuing use of commutative 
algebra and algebraic geometry, as well as operator algebra and C * -algebra tech- 
niques. 



4. Operator theory in spaces with an indefinite metric 

Recent achievements. Indefinite inner products have arisen in several natural ways 
over the last couple of decades: J-inner matrix functions arising in the parametriza- 
tion of all solutions of a matricial Nevanlinna-Pick interpolation or of a Nehari 
problem, solution of the inverse problem of characterizing which J-inner func- 
tions parametrize all solutions of an interpolation problem, the Adamjan-Arov- 
Krein generalization of Nehari’s theorem to handle the case of approximation of 
an L°°- function by a function in the class H £° (bounded on the unit circle and 
meromorphic on the unit disk with at most k, poles), extensions of the concept of 
reproducing kernel Hilbert space to the notion of reproducing kernel Pontryagin 
space. A fundamental result now relatively old is the modified spectral theorem 
for an operator selfadjoint in a Pontryagin space. Recent achievements in this 
area represented at IWOTA2002 include: interpolation problems involving Pick 
matrices having some number k of negative eigenvalues, spectral theory for spe- 
cial cases of matrices normal in an indefinite inner product, spectral properties of 
polynomials orthogonal in an indefinite inner product, spectral theory of canonical 
differential equations involving an indefinite inner product, Darlington representa- 
tions in an indefinite-inner-product setting, an analysis of the Schur algorithm for 
a function in the generalized Schur class S K , new approaches to and understand- 
ing of the Krein-Langer factorization for functions in the class S K , and continuing 
exploitation of Krem-space operator theory techniques for the analysis of assorted 
generalized boundary- value problems for differential operators coming up in engi- 
neering and mathematical-physics contexts. 

Future directions. There are many concrete problems in differential equations or 
mathematical physics where the theory of operators in indefinite metric spaces 
is relevant. Examples include the problem of the motion of an ideal or viscous 
liquid around a solid body, the theory of linear systems having Pontryagin state 
spaces, and direct and inverse problems for indefinite cases of canonical differential 
equations. We expect to see continuing interest in and work on such applications 
of Krein-space operator theory. On the more theoretical side, development of a 
structure theory for various new classes of operators in indefinite metric spaces, 
even in finite-dimensional spaces, remains a challenge. 
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5. Operator theory in harmonic analysis 

Recent achievements. Recent achievements include the introduction of operator- 
theoretic techniques, in particular, the optimality principle for the principle of 
dynamic programming in control theory, for new treatments and new solutions of 
problems in harmonic analysis (e.g., estimates for Carleson measures and Hunt- 
Muckenhoupt-Wheeden weights, including new results for the case of matrix- or 
operator- valued weights). 

Recently X. Tolsa solved the Painleve and Vitushkin problems on analytic 
capacity. The major contribution of Vitushkin was the reduction of problems in ra- 
tional approximation to analytic capacity estimates, which then made verification 
of subadditivity conditions crucial. Now it is possible to give an equivalent, more 
operator-theoretic formulation of analytic capacity, thereby making subadditivity 
easier to prove. 

Future directions. The full operator-theoretic reformulations and implications of 
the recent solution of Vitushkin’s problem remain to be explored. There should 
follow more explicit results in rational approximation and harmonic analysis in 
which operator theory will play a prominent role. 



6. Factorization of matrix functions 

Recent achievements. Factorization of matrix functions was developed in the first 
half of the last century as a tool for studying and solving systems of singular 
integral equations or, equivalently, certain boundary value problems for analytic 
functions. In the 70s and 80s its applications to inverse scattering theory were 
found and intensively developed. Recently researchers have discovered several more 
rather unexpected applications, including the theory of random matrices, asymp- 
totics of orthogonal polynomials, completion and extension problems for (almost) 
periodic functions of several variables. As a result, some new approaches were de- 
veloped in the factorization theory per se, such as the steepest descent method in 
the asymptotical analysis of parametrized factorization problems. The cornerstones 
of the intrinsic development of the factorization theory also include the complete 
treatment of the piecewise continuous case on Carleson curves with Muckenhoupt 
weights, and the factorization theory of semi-almost periodic matrix functions. 

Future directions. Several classical problems of factorization theory still have ap- 
peal. We mention here the problem of explicit factorization of new classes of matrix 
functions, and we expect further progress to be achieved here with more systematic 
use of special functions and the Riemann boundary problem on Riemann surfaces. 
Computation of partial indices and, in particular, finding new verifiable conditions 
for their stability, remains an important open problem. Substantial progress in this 
direction was achieved lately using the current ideas of numerical methods, but 
more remains to be done. 
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New problems, arising from the needs of the latest developments (both theo- 
retical and applied), include the systematic treatment of matrix functions defined 
on curves with intersections and the asymptotic behavior of the factorization in 
case of parameterized families of matrices more general than piecewise continuous 
(in particular, having discontinuities of semi-almost periodic type). Though the 
properties of the almost periodic factorization, when it exists, are rather well un- 
derstood by now, the intrinsic nature of its existence remains a mystery and calls 
for further research. Work on explicit factorization of concrete almost periodic ma- 
trix functions, arising from convolution type equations on finite intervals, is also in 
order. The recently discovered connections with the Corona problem might prove 
to be useful in this regard. 

Yet another natural direction of development is the factorization theory of 
matrix functions defined on abstract ordered abelian groups. Such a theory, when 
constructed, would embrace the classical results on Wiener-Hopf factorization of 
continuous matrix functions on the unit circle, on the one hand, and (already 
mentioned) more recent almost periodic factorization, on the other. One should 
anticipate completely new phenomena here, and a need to develop radically new 
ideas and approaches; indeed, even the almost periodic factorization is far from 
being analogous to the classical case of functions defined on the unit circle. 



7. Integration of operator theory with wavelets 

Recent achievements. While the theory of wavelets and its applications to data 
compression and image reconstruction has been undergoing intense development 
in the past couple of decades, the understanding of its connections with opera- 
tor theory (and operator algebra) ideas is more recent (in the work of Larson, 
Jorgensen and their collaborators), and the interaction (between wavelets and op- 
erator theory) goes both ways. Here are a few samples of the connections: 

The theory of wavelets has led to a generalization of Halmos’s idea of a 
wandering vector for a (unitary or isometric) shift operator on a Hilbert space to 
the notion of wandering vector for a unitary system (with generators satisfying 
some sort of intertwining condition rather than commuting) . This operator-theory 
point of view has then given further insight and new results in the theory of 
wavelets. 

A fundamental operator for many wavelet considerations has turned out to 
be the so-called transfer operator. Indeed, this operator determines the regularity 
properties of wavelets, e.g., the strength of the approximations, and additional 
fine structure of direct practical significance. The spectral theory of this opera- 
tor, including application of an infinite-dimensional extension of Perron- Frobenius 
theory, has led to better understanding of wavelet constructions. The transfer op- 
erator has applications in statistics, mathematical physics, and dynamical systems. 
Within pure mathematics, it is studied in analytic number theory and harmonic 
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analysis in connection with zeta functions and trace formulas. The spectral the- 
ory of the operator varies with the context of the function spaces: these are the 
regularity spaces which, at the same time, play a crucial role in questions around 
discretization of PDE problems (e.g., finite element methods). The operator the- 
oretic aspects of wavelet analysis are covered strongly at a number of centers, in 
the US, Europe, and Asia. The wavelet algorithm is now known to have a very fast 
polynomial quantum version, while the classical one is not known to be polynomial. 
This is an area with a lot of interest in the physics community. 

Future directions. We expect to see more interaction between the wavelet and op- 
erator theory communities. Wavelet analysis has led to a number of problems in 
areas which have traditionally been in the domain of the operator theory commu- 
nity, e.g., spectral theory of matrix polynomials, factorization of unitary-operator- 
valued functions, Fredholm index and related homotopy questions, spectral and 
index theory for concrete operators on a function space. To this point, researchers 
in the wavelet community have generally come up with their own ad hoc solutions 
rather than using the systematic theories developed in the operator-theory commu- 
nity. On the other hand, the problems in wavelet theory are often slightly different 
than anything which has already been done in the operator-theory community. For 
example, the transfer operator mentioned above does not fit into established trends 
in operator theory - the closest, but not exact, fit is perhaps the Cowen-Douglas 
class of operators. These various connections give a fruitful area for continuing 
cross- fertilization. 



8. Structured matrices 

Recent achievements. Structured matrices, such as Toeplitz, Hankel, Vander- 
monde, Cauchy, Pick, Bezout, separately and in combinations, such as Toeplitz- 
plus-Hankel, or more generally, matrices with “displacement structure” , have been 
extensively studied for a long time, from diverse and seemingly unrelated points 
of view, in mathematics, computer science, and engineering. This includes also the 
point of view of operator theory: for example, Toeplitz matrices have been stud- 
ied using methods of reproducing kernel Hilbert space and lifting of commutants. 
A relatively recent achievement is the development of fast solution algorithms 
for Toeplitz, Vandermonde and Cauchy matrices, solutions of Lyapunov/ Riccati 
equations, types of Pick matrices. One approach to these fast solution algorithms 
is based on the theory of rational matrix interpolation. In particular, there now 
exist superfast algorithms for solving passive tangential interpolation problems 
of Nevanlinna-Pick type, and a new algorithm for list decoding of Reed-Solomon 
codes. New ideas entering the mainstream are the notion of pseudospectrum and 
applications of spectral theory for various classes of infinite-dimensional operators, 
and of random matrices to understanding conditioning and asymptotics issues for 
various Galerkin approximations. 
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Future directions. There remain matrices for which fast algorithms should be de- 
veloped: circulant matrices, multi-dimensional block Toeplitz (existing algorithms 
do not use the complete special structure). We expect to see more interaction be- 
tween numerical analysis, the theory of random matrices, and the spectral theory 
of infinite-dimensional operators. 



9. Linear matrix inequalities 

Recent achievements. The recent development of efficient interior-point algorithms 
for convex optimization problems involving linear matrix inequalities (LMIs) has 
spurred research in a wide variety of application fields, including control system 
analysis and synthesis, combinatorial optimization, circuit design, structural opti- 
mization, experiment design, and geometric problems involving ellipsoidal bound- 
ing and approximation. The interior point methods for convex optimization were 
reviewed in a plenary lecture at IWOTA2002. 

Future directions. Linear system problems convert readily to systems of matrix 
inequalities; then the trick (often impossible) is to convert them to “nice” matrix 
inequalities. There are only a few tricks known, the main one being to use “Schur 
complements” . Despite the attention of hundreds of engineers, essentially nothing 
general or systematic is known. Matrix inequalities fall directly within the exper- 
tise of the operator-theory community; recently engineering-type inequalities are 
getting scrutiny in this community. A major issue central to numerical solution is 
which problems are convex matrix inequalities and which convert to linear matrix 
inequalities. There are beginning to be results on such questions and the pace is 
quickening. 



10. Distributed-parameter system theory 

Recent achievements. Axiomatization of the integral form of an input-state-output 
system (involving bounded operators between input, state and output signals over 
each time interval rather than between instantaneous input, state and output) in 
the notion of “well-posed system” has led to formulation of a useful rigged setup for 
an operator colligation with unbounded infinitesimal input, state dynamics, and 
output operators. Such colligations (including extensions corresponding to sys- 
tems which are not well-posed) have led to new realization theorems for functions 
which are positive-real, Herglotz-Nevanlinna, or bounded-real over a half-plane. 
This foundation in turn is leading to better understanding of state-space methods 
for feedback stabilization, J-spectral factorization, and optimal and iif^-control 
for distributed parameter systems. There have been also interesting connections 
between function theory (e.g., the corona problem) and well-posedness questions 
for a given dynamics and unbounded input or output operator. 
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Future directions. Application of the theory to new concrete examples and appli- 
cations, including such situations as quantum systems and linearized fluid flow, is 
an important future direction for the theory. Also integration of the more general 
behavior-theoretic/implicit-system viewpoint can be expected. 



11. Direct and inverse scattering and spectral problems 

Recent achievements. Riemann-Hilbert problems/ Wiener-Hopf factorization: In- 
verse scattering problems arise in wave propagation and quantum mechanics: one 
seeks to recover the potential function or the differential operator from knowledge 
of the scattering data or spectral measure. The direct problem is to compute the 
scattering function and spectral measure from the operator. Well-known inverse 
spectral methods are via solving Riemann-Hilbert problems and related integral 
equations and Wiener-Hopf factorization problems, or by making use of the the- 
ory of Riemann surfaces and theta functions. A more recent method borrows the 
so-called state-space method from system theory: one assumes that the scattering 
function is rational and then solves the inverse scattering problem in terms of the 
matrices arising from the realization of the reflection coefficient as the transfer 
function of a linear system. The connection with interpolation problems is ex- 
ploited to obtain solutions of bitangential inverse scattering and spectral problems 
for canonical integral and differential systems. Formulas are obtained via reproduc- 
ing kernel Hilbert space methods with kernels based on strongly regular J-inner 
matrix- valued functions. 

Nonlinear PDEs: Operator theory helps to analyze various linear and non- 
linear differential equations describing important physical phenomena. One im- 
portant application is the inverse scattering transform; specifically, solutions of 
certain nonlinear partial differential equations (nonlinear evolution equations) are 
obtained by solving the inverse scattering problems of related linear ordinary differ- 
ential equations. Operator-theory techniques to solve inverse scattering problems 
thus play an important role in solving various nonlinear partial differential equa- 
tions. Some special solutions of such nonlinear partial differential equations, known 
as soliton solutions, have led to advances in communications via fiber optics. The 
analysis of the inverse scattering transformation leads to questions in the spectral 
theory of nonselfadjoint systems of differential equations, like the Zakharov- Shabat 
system or the Manakov system. Research in this area has been of interest for both 
mathematicians and engineers. 

Ergodic and almost ergodic operators: Functional models of operators on Rie- 
mann surfaces which were introduced in operator theory are now being used to de- 
velop spectral and scattering theories of ergodic and almost periodic difference/dif- 
ferential operators. 

Future directions. Perturbations of selfadjoint (and normal) operators which are 
not, say, selfadjoint or even accretive: Applications include spectral problems aris- 
ing in the theory of pulse propagation in nonlinear optical fibers. 
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Selfadjoint perturbations in Schatten-von Neumann ideals larger than the 
trace ideal, and Schrodinger operators with long range non-smooth potentials: The 
two issues mentioned here are connected and actually related to an old question 
of M.G. Krem raised at the Moscow Mathematical Congress (1966): find spec- 
tral interpretations of M.M. Djrbashian’s factorization of meromorphic functions 
from generalized Nevanlinna classes. The connection between the two is as follows. 
Schrodinger operators with long-range potentials can be viewed as non-trace class 
perturbations of free Hamiltonians. Consequently, many analytic functions (like 
the perturbation determinant) appearing in the trace class perturbation theory 
need suitable modifications. The problem is that Nevanlinna functions do not be- 
long to this setting. There are strong indications to conjecture that the relevant 
functions are from the so-called Djrbashian class of analytic functions of gener- 
alized bounded type but it appears to be a very difficult issue. Thus, important 
developments in this direction should also address the Krem problem. 

A proper solution of questions arising in the spectral theory of almost periodic 
difference/differential operators alluded to above requires a real breakthrough in 
the subject, namely, functional models on Riemann surfaces of infinite genus going 
beyond the Widom class. 
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Inverse Scattering Transform, KdV, 
and Solitons 

Tuncay Aktosun 



Abstract. In this review paper, the Korteweg-de Vries equation (KdV) is con- 
sidered, and it is derived by using the Lax method and the AKNS method. An 
outline of the inverse scattering problem and of its solution is presented for 
the associated Schrodinger equation on the line. The inverse scattering trans- 
form is described in order to solve the initial- value problem for the KdV, and 
the time evolution of the corresponding scattering data is obtained. Soliton 
solutions to the KdV are derived in several ways. 



1. Introduction 



The Korteweg-de Vries equation (KdV, for short) is used to model propagation 
of water waves in long, narrow, and shallow canals. It was first formulated [1] in 
1895 by the Dutch mathematicians Diederik Johannes Korteweg and Gustav de 
Vries. Korteweg was a well known mathematician of his time, and de Vries wrote 
a doctoral thesis on the subject under Korteweg. 

After some scaling, it is customary to write the KdV in the form 



du du d 3 u 
dt ~ + !h? =0. 



ieR, t> o, 



(l.i) 



where —u(x,t) corresponds to the vertical displacement of the water from the 
equilibrium at the location x at time t. Replacing u by — u amounts to replacing 
—6 by +6 in (1.1). Also, by scaling x, £, and u, i.e., by multiplying them with 
some positive constants, it is possible to change the constants in front of each of 
the three terms on the left-hand side of (1.1) at will. 

Note that the KdV is a nonlinear partial differential equation (PDE, for 
short) due to the presence of the uu x term, where we use a subscript to denote 
the partial derivative. The u xxx term makes it dispersive , i.e., in general an initial 
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wave u(x,0) will broaden in space as time progresses. In addition to its solutions 
showing behavior of nonlinearity and dispersiveness, the KdV possesses certain 
special solutions, known as solitary wave solutions, which would not be expected 
from a nonlinear and dispersive PDE. A single solitary wave solution to a PDE 
retains its shape in time and moves either to the left or right. It turns out that 
such a solution to (1.1) only moves to the right, and any general solitary wave 
solution asymptotically resembles a train of single solitary wave solutions. 

It was a Scottish engineer named John Scott Russell who first observed a 
solitary water wave. This happened in 1834 on the Edinburgh- to- Glasgow canal, 
some 60 years before the formulation of the KdV as a model for water waves. 
Russell reported [2] his observation to the British Association for the Advancement 
of Science in September 1844. The details of Russell’s observation, including the 
rivalry between him and George Airy who did not believe in the existence of 
solitary water waves, can be found in [3-6]. 

As outlined in the famous preprint [7], which was never published as a jour- 
nal article, Enrico Fermi in his summer visits to Los Alamos, jointly with J. Pasta 
and S. Ulam, studied a one-dimensional (1-D, for short) dynamical system of 64 
particles with forces between neighbors containing nonlinear terms. These compu- 
tational studies were carried out on the Los Alamos computer named Maniac I. The 
primary aim of the study was to determine the rate of approach to the equiparti- 
tion of energy among various degrees of freedom. Contrary to expectations, Fermi, 
Pasta, and Ulam observed very little, if any, tendency towards the equipartition 
of energy, but instead the almost ongoing recurrence to the initial state, which 
was puzzling. After Fermi’s death in November 1954, Pasta and Ulam completed 
their last few computational examples. Their preprint appears in Fermi’s Collected 
Papers [8] and is also available on the internet [9] . 

The importance of the KdV arose in 1965, when Zabusky and Kruskal [10] 
were able to explain the Fermi-Pasta-Ulam puzzle in terms of solitary- wave solu- 
tions to the KdV. In their analysis of numerical solutions to the KdV, Zabusky 
and Kruskal observed solitary-wave pulses, named such pulses solitons because of 
their particle-like behavior, and noted that such pulses interact with each other 
nonlinear ly but come out of their interaction virtually unaffected in size or shape. 
Such unusual nonlinear interactions among soliton solutions to the KdV created 
a lot of excitement, but at that time no one knew how to solve such a nonlinear 
PDE, except numerically. 

In their celebrated paper [11] of 1967, Gardner, Greene, Kruskal, and Miura 
presented a method, now known as the inverse scattering transform, to solve 
the initial- value problem for the KdV, assuming that the initial value u(x,0) ap- 
proaches a constant sufficiently rapidly as x —> ±oo. There is no loss of generality 
in choosing that constant as zero. They showed that u(x,t) can be obtained from 
u(x , 0) with the help of the solution to the inverse scattering problem for the 1-D 
Schrodinger equation with the time-evolved scattering data. They also explained 
that soliton solutions to the KdV correspond to the case of zero reflection coeffi- 
cient in the scattering data. They observed from various numerical studies of the 
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KdV that, for large t, u(x,t) in general consists of a finite train of solitons trav- 
eling in the positive x direction and an oscillatory train spreading in the opposite 
direction. 

In our paper, we present an elementary review of the inverse scattering trans- 
form for the KdV. We consider the time-evolved Schrodinger equation, where V ( x ) 
in (2.1) is replaced by u(x,t); namely, we deal with 

- ^^ 2 ^ + x \ t) — u ( x > t) x 'it)i x G R, (1.2) 

where t > 0 is a parameter that is usually interpreted as time. Thus, we view V (x) 
as the initial value u(x, 0) of the potential u(x, t) and look at ^>(fc, x; t) as the time 
evolution of ^(fc,x) of (2.1) from the initial time t = 0. The scattering coefficients 
T(k\t), R(h,t), and L{k\t ) associated with (1.2) are viewed as evolving from the 
corresponding coefficients T(fc), R(k), and L(k) of (2.1), respectively, from t = 0. 
Thus, our notation is such that 

V(x) = u(x, 0), ^>(fc,x) = ^(fc,x;0), 

T(k) = T(fc; 0), R(k) = R(k ; 0), L(k) = L(fe; 0). 

Our review paper is organized as follows. In Section 2 we consider scattering 
solutions and bound state solutions to the Schrodinger equation (2.1) and intro- 
duce the scattering coefficients, bound- state norming constants , and dependency 
constants corresponding to a potential in the so-called Faddeev class. In Section 
3 we present an outline of the inverse scattering problem for (2.1) and review 
some solution methods based on solving an associated Riemann- Hilbert problem. 
In Section 4 we consider the time evolution of the scattering coefficients, bound- 
state norming constants, and dependency constants when the potential evolves 
from u(x,0) to u(x,t ). We also introduce the Lax pair associated with the KdV 
and derive the KdV by using the Lax method. In Section 5 we study the AKNS 
method and derive the KdV via that method. In Section 6 we present some of the 
methods to solve the initial value problem for the KdV. In Section 7 we concentrate 
on soliton solutions to the KdV and obtain various representations of the AT-soliton 
solution. Finally, in Section 8 we provide certain remarks on the Backlund trans- 
formation, the conserved quantities, and some other aspects related to the KdV. 



2. Schrodinger equation and the scattering data 

Consider the Schrodinger equation 

~^dx^ ' = V(x)'ip(k,x), x G R, (2.1) 

where V is real- valued and belongs to L\{ R). Here, L^(R) denotes the class of 
measurable potentials such that dx (1 + |x| n ) | V(x)\ is finite. The class of real- 
valued potentials in L}(R) is sometimes called the Faddeev class, after Ludwig 
Faddeev’s analysis [12] of the inverse scattering problem for (2.1) within that 
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class of potentials. In appropriate units, (2.1) describes the quantum mechanical 
behavior of a particle of total energy k 2 under the influence of the potential V. 
The inverse scattering problem for (2.1) consists of the determination of V from 
an appropriate set of scattering data. 

There are two types of solutions to (2.1). The scattering solutions consist of 
linear combinations of e lkx and e~ lkx as x — > ±oo, and they occur for k G R\ {0}. 
Real k values correspond to positive energies of the particle, and a particle of 
positive energy can be visualized as being capable of escaping to x = d=oo as a 
result of scattering from V. Heuristically, since V(x) vanishes at x = ±oo, the 
particle will still have some kinetic energy at infinity and hence is allowed to be 
at infinity. On the other hand, a bound state of (2.1) is a solution that belongs 
to L 2 (R) in the x variable. It turns out that, when V belongs to the Faddeev 
class, the bound-state solutions to (2.1) decay exponentially as x — * ±oo, and 
they can occur only at certain A;- values on the imaginary axis in C + . We use C + 
to denote the upper-half complex plane and put C+ := C + UR. Each bound state 
corresponds to a negative total energy of the particle, and as a result the particle 
is bound by the potential and does not have sufficient kinetic energy to escape to 
infinity. We will use N to denote the number of bound states, which is known to 
be finite when V is in the Faddeev class, and suppose that the bound states occur 
at k = %Kj with the ordering 0 < K\ < • • • < kn- 

Among the scattering solutions to (2.1) are the Jost solution from the left, 
/i, and the Jost solution from the right, / r , satisfying the respective boundary 
conditions 



e~ ikx fi(k,x) = l + o(l), e~ ikx f{(k,x) =tA + o(l), x -> +oo, (2.2) 

e ikx fr(k,x) — 1 + o(l), e ikx f'(k,x) = -ik + o(l), x -» -oo, (2.3) 

where the prime is used for the derivative with respect to the spatial coordinate 
x. From the spatial asymptotics 

p ikx T(k\p~^ x 

*M-W) + m- +all) - <2 ' 4) 

p —ikx R(lc\ p^kx 

mx )= m + T ^)~ +0{11 (2 ' 5) 

we obtain the scattering coefficients, namely, the transmission coefficient T, and 
the reflection coefficients L and R , from the left and right, respectively. It is also 
possible to express the scattering coefficients in terms of certain Wronskians [12- 
16] involving f\ and / r . We have 



T(k) = 



2 ik 

[Mk,x);fi(k,x)]’ 



T(!A [fi(k,x); f r (—k,x)] [fi(-k,x);f v (k,x)\ 

[fr(k,x);fi(k,x)\ ’ ^ [fr(k,x);fi(k,x)] 

where the Wronskian is defined as [F; G] := FG' — F'G. 



( 2 . 6 ) 
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It is known [12-16] that, for each fixed x G R, the Jost solutions /i(-,x) and 
/ r (*,x) have analytic extensions in k to C + . Moreover, 

fi(-k*,x) = fi(k,x)*, f r (-k*,x) = f r (k,x)*, k e C+, 

T(-k) = T(k)*, R(-k) = R(k)*, L(-k) = L(k)*, ke R, 

where the asterisk denotes complex conjugation. We also have 

R(k) T(k)* = -L(k)* T(k), ke R, (2.7) 

\T(k)\ 2 + \L(k)\ 2 = 1 - \T(k)\ 2 + \R{k)\ 2 , (2.8) 

Thus, the scattering coefficients cannot exceed one in absolute value for real k. 
Furthermore, T(k) ^ 0 if k G R \ {0}, and hence the reflection coefficients are 
strictly less than one in absolute value when k G R \ {0}. In general, R and L 
are defined only for real k values, but T has a meromorphic extension to C + . For 
large k one has 

T(k) = 1 + 0(l/k), k — » oo in C+, 

R(k) = o(l/k), L(k) = o(l/fc), k ± oo. 

Each bound state corresponds to a pole of T in C + and vice versa. It is 
known that the bound states are simple, i.e., at each k = %Kj there exists only 
one linearly independent solution to (2.1) belonging to L 2 (R). The bound-state 
norming constants c\j and c r j are defined as 

1 - 1/2 





" roo 


-1/2 


z -00 


c\j 


/ dx f\(iKj,x) 2 

.J —OO 


, c r j := 


/ dxf r (iKj,xY 

— OO 



and they are related to each other via the residues of T as 

c 2 - 

Res (T, iKj) = i c 2 j 7 j = i — , 
where 7 j is the dependency constant given by 



7 j : ~~ 






c 2 - 
7 j 


(2.9) 




(2.10) 



The sign of 7 j is the same as that of (— 1)^ i and hence 



-rj 



-(-1) 



N-j 






The normalized bound-state solution <fij(x ) at k = %Kj is defined as 
<fij(x) := Cij f\{iKj,x) = (-1 ) N ~ J c rj f r (iKj,x). 
The scattering matrix associated with (2.1) is given by 

'T(fc) R(k) 

L(k) T(k) ’ 



S (k) := 



ke R, 
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and it can be constructed in terms of the bound-state energies and either one of 
the reflection coefficients R and L. Given R(k) for k G R and the bound-state 
poles k = ik j, one can construct T as 



T(k) = 




log(l - I R(s)\ 2 ) 

s — k - i0+ 



keC+, (2.11) 



and use (2.7) to construct L(k) for k e R. Similarly, given L(k) for k G R and the 
bound-state poles k = iKj, one can construct T [cf. (2.8) and (2.11)] as 



T(k) = 




log(l ~ l£(s)l 2 ) 

s - k — i0+ 



k £ C+, 



and obtain R(k) for k G R via (2.7). 



3. Inverse scattering problem 

When there are no bound states, either one of the reflection coefficients R and L 
uniquely determines the corresponding potential in the Faddeev class. However, 
when there are bound states, for the unique determination of V, in addition to 
one reflection coefficient and the bound-state energies, one must also specify a 
bound-state norming constant or, equivalently, the dependency constant for each 
bound state. To recover V uniquely, as our scattering data we may use either the 
left scattering data {R, {ftj} , {c\j}} or the right scattering data {L, {^}, {c rj }}; 
these two are equivalent to each other, and each is also equivalent to {S, {7^}}. 

A characterization for a specific class of potentials consists of specifying some 
necessary and sufficient conditions on the scattering data which guarantee that 
there exists a corresponding unique potential in that class. Such conditions are 
usually obtained by using the Faddeev- Marchenko method [12-19]; this method is 
also known as the Marchenko method, and sometimes in the literature it is referred 
to as the G el Jand- Levitan- Marchenko method even though this is a misnomer [20]. 
The characterization conditions can be stated for the left scattering data, for the 
right scattering data, or for the combination of both. For a characterization in the 
class of real- valued potentials belonging to R), the reader is referred to [14]. 
Various characterizations in the Faddeev class can be found in [13,17-19]. 

Since k appears as k 2 in (2.1), the functions f\(—k,x ) and f r (—k,x ) are also 
solutions to (2.1) and they can be expressed as linear combinations of the Jost 
solutions f\(k,x) and f r (k,x) as 

f\(-k, x) = T{k) f r (k, x) - R{k) f\(k , x), k G R, 

/ r (-fc, x) = T(k) /i(fc, x) - L(k) f T (k, x), ke R, 

or equivalently as 

m\(—k, x) = T(k ) m r {k , x) — R(k) e 2lkx m\(k , x), 



k e R, (3.1) 
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m r (—k,x)=T(k)m\(k 1 x)—L(k)e 2lkx m r (k 1 x), k G R, (3.2) 
where m\ and ra r are the Faddeev functions defined as 

m\(k,x) := e~ lkx f\{k,x), m r (k,x) := e lkx f r (k,x). (3.3) 

Each of (3.1) and (3.2) can be viewed as a Riemann-Hilbert problem [6,15], where, 
knowing the scattering coefficients for k G R, the aim is to construct m\ and m r 
such that, for each xgR, m\(-,x) and m r (*,x) are analytic in C + , continuous in 
C+, and behave like 1 + 0(1/ k) as k — > oo in C + . Once m\(k,x) or m r (fc,x) is 
constructed, the potential can be obtained with the help of (2.1) and (3.3), namely, 
by using 



V(x) 



m['(k,x) m((/c,x) 

rri\(k,x) 1 m\(k,x) J 



V(x) 



m"(k,x) _ 2i] z m T(k,x) 
m r (k,x) m r (k,x)' 



(3.4) 



where the right-hand sides can be evaluated at any particular value of k G C+. 

Alternatively, the potential can be constructed by the Faddeev- Marchenko 
method; namely, V can be obtained from the left scattering data {R, {/%•}, {c\j}} 
by solving the left Marchenko integral equation or from the right scattering data 
{L, {ftj}, {c r:/ }} by solving the right Marchenko integral equation. 

The left Marchenko equation using the left scattering data as the input is 
given by 

f‘00 

B\(x,y) + Q](2x + y) + / dz Cl\(2x + y + z) B\(x, z) = 0, y > 0, (3.5) 

Jo 



where 



fli (y) :=T r dkR(k), 

J — oo 



2 --Kjy 



One can obtain (3.5) from (3.1) via a Fourier transformation. Once (3.5) is solved 
and B\(x,y) is obtained, the potential is recovered as 



V(x) = -2 



dBi(x,0+) 



and the Faddeev function from the left is constructed as 

/*OC 

rri](k,x) = l+ dy Bi(x,y) e tky . (3.7) 

Jo 

Similarly, via a Fourier transformation on (3.2), using the right scattering 
data as the input one obtains the right Marchenko equation 

poo 

B r (x,y) + fi r (— 2x + y) + / dz O r (— 2x 4- y + z) B r (x, z) — 0, y > 0, (3.8) 

Jo 



where 



-j poo N 

:= ^ J_ dk L(k) e iky + £ c 2 rj e 
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Once (3.8) is solved, the potential is recovered by using 

ax 

and the Faddeev function from the right is constructed as 



POO 

m T (k,x) = l+ dyB r (x,y)e lky . (3.10) 

Jo 

When the characterization conditions on the scattering data corresponding to po- 
tentials in the Faddeev class are satisfied, both the left and right Marchenko equa- 
tions are uniquely solvable, and the right-hand sides of (3.6) and (3.9) are equal 
to each other and belong to L\(R). Thus, V can be obtained from either (3.6) 
or (3.9). There are various other methods to recover V from an appropriate set 
of scattering data. We refer the reader to [16,19] for a review of some of those 
methods. 



4. Lax method and evolution of the scattering data 

Soon after Gardner, Green, Kruskal, and Miura demonstrated [11] that the initial 
value problem for the KdV can be solved by the inverse scattering transform, Peter 
Lax gave [21] a criterion to show that the KdV can be viewed as a compatibility 
condition related to the time evolution of solutions to (1.2). Since Lax’s criterion 
is applicable to other nonlinear PDEs solvable by an inverse scattering transform 
(e.g., the initial- value problem for the nonlinear Schrodinger equation can be solved 
[4,6,22] via the inverse scattering transform for the Zakharov-Shabat system), we 
first outline the general idea behind the Lax method and next demonstrate its 
application on the KdV. 

Given a linear operator C with C'tp = A^, we are interested in finding another 
operator A (the operators A and C are said to form a Lax pair) such that: 

(i) The spectral parameter A does not change in time. 

(ii) The quantity ^ — Ai/j must remain a solution to C'i/j = Xip. 

(iii) The quantity C t + CA — AC must be a multiplication operator. 

As the following argument shows, for compatibility, we are forced to have 

C t + CA-AC = 0, (4.1) 



which is interpreted as an integrable PDE and in general is nonlinear. From con- 
dition (ii) above we see that 

C ($ t - A&) = A ('ipt ~ A^) , 

or equivalently, 

C'lpt — CA'ip = X'lpt — A (X'tp) 

= d t (X'tp) - AC'ip 
= d t (Cijj) - AC'il) 

= C t % b + Ctyt ~ AC^jj, 



(4.2) 
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where we have used C ^ = A ^ and X t = 0. After canceling the Cipt terms in (4.2), 
we get 

(C t + CA - AC) ^ = 0. (4.3) 

Because of condition (iii) listed above, from (4.3) we obtain the compatibility 
condition (4.1). 

Let us write the Schrodinger equation (1.2) as 

Cif) = A-0, C := —d x + u(x, t), (4.4) 



and try 



A = ad? + 0dl + £d x + r?, 



(4.5) 



where the coefficients a, /?, £, and rj may depend on x and £, but not on A. Note 
that Ct = ut and A is the same as k 2 . Using (4.4) and (4.5) we can write (4.1) 
explicitly as 

(R 5 + (R 4 + ()^ + (R 2 + ()a x + ()=0, (4.6) 

where each ( ) in (4.6) denotes the appropriate coefficient. The coefficient of 
automatically vanishes. Setting the coefficients of d 3 x equal to zero for j = 4, 3, 2, 1, 
we obtain 



a = ci, /3 = c 2 , £ = c 3 --ciu, rj = C4 - -ciu x - c 2 u, 

where c i, C2, c 3 , and c\ are arbitrary constants. Using c\ = —4 and c 3 = 0 in the 
last term on the left-hand side of (4.6) and setting that last term equal to zero, we 
obtain the KdV equation given in (1.1). Further, by using C2 = C4 = 0, we obtain 
the operator A associated with C as 



A = -4<9f + 6ud x + 3 u. 



(4.7) 



Let us remark that condition (ii) stated above allows us to determine the 
time evolution of any solution to the Schrodinger equation (1.2) as the initial 
potential u(ar, 0) evolves to u(x,t). For example, let us find the time evolution of 
/i(A:,x;t), the Jost solution from the left. By using (2.2) and (2.5) with T(k) and 
R(k) replaced by T(k;t) and R(k;t ), respectively, we see that 

e~ ikx fi{k,x- : t) = l+o(l), e- ikx f[(k,x-,t) =ik + o(l), x -> +oo, (4.8) 



f r {k,x;t) 



Q — ikx 



+ 



R{k\ t ) e 



ikx 



+ 0 ( 1 ), 



T(k;t) T(k; t) 

From condition (ii) of the Lax method and (4.7) we obtain 



+oo. 



(4.9) 



d t fi(k,x;t) - (-4 dl+ 6ud x + Su x ) fi(k,x;t) 

= p(k, t) fi(k, x; t ) + q(k, t) f r (k, x; t), 



(4.10) 



where we have used the fact that the quantity d t fi — Afi remains a solution to 
(1.2) and hence can be expressed as a linear combination of the two linearly inde- 
pendent Jost solutions fi(k, x; t ) and f r (k, x; t) with coefficients p(k, t) and q(k, t), 
respectively. For each fixed t, assuming that u(x,t) = o(l) and u x (x,t) = o(l) as 
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x — > +oo, the coefficients p(k,t ) and q(k,t) can be evaluated by letting x — > +oo 
in (4.10). Using (4.8) and (4.9) in (4.10), we get 



d t e ikx + Adle ikx = p(k, t) e lkx + q(k, t) 



4 —ikx . ikx 

T(k;t) T{k;t) 



(4.11) 



Prom (4.11), by comparing the coefficients of e lkx and e lkx on both sides, we 
obtain 



q(k,t) — 0, p(k,t) = —Aik 3 . (4-12) 

Thus, the time evolution of /i(fc, x ; t) is determined by the linear third-order PDE 
d t fi(k, x; t) - Af](k, x; t) - -Aik? f\{k, x; t ). (4-13) 



Similarly, letting x — > — oo in (4.10), with the help of (2.3), (2.4), and (4.12) we 
obtain 



ft 



1 

T(M) 



e ikx + 



L(k; t) ik x 

T(k; t) 



= [- 4^ 3 



4iA: 3 



^ ikx . —ikx 

T(k-t) T(M) 



(4.14) 



where we have also used that u(x, t ) = o(l) and u x (x, t) = o(l) as x — > — oo. From 
(4.14), comparing the coefficients of e lkx and e~ lkx on both sides, we obtain 



ftT(ife; t) = 0, d t L{k- 1) = —8ik 3 L(k; t), 



and hence 

T(fc; t) = T(k\ 0) = T(/c), L(k; t) = L(fc; 0) e~ 3ikH = L(k) e~ 3ikH . (4.15) 

Thus, the transmission coefficient remains unchanged and the reflection coefficient 
from the left undergoes a simple phase change as t progresses. 

Proceeding in a similar manner, we can obtain the time evolution of the Jost 
solution f T (k,x;t) and the right reflection coefficient R(k;t). With A as in (4.7), 
we get 

dtfi(k, x; t) - Af r (k, x; t ) = 4 ik 3 f r (k, x\ t), (4.16) 

R(k; t) = R(k ; 0) e 8lfc3 * = R(k) e* ikH . (4.17) 

In order to evaluate the time evolution of the dependency constants 7 j(t), 
we can substitute 7 j(t)f T (iKj,x]t) for f\(iKj,x',t) [cf. (2.10)] and evaluate (4.13) 
at k = iKj. We get 

f r (iKj,x; t ) daj{t) + 7 j(t) dtfr(iKj,x; t)-'y j (t)Af T (iK j ,x; t) 

On the other hand, from (4.16) at k — iKj, we obtain 

7j(t)dtf T (iKj,x;t) -'y j (t)Af r (iK j ,x;t) = 4«; 3 7 J (f)/ r (m J -,x;f). 

Subtracting (4.19) from (4.18) we conclude that dtjj(t) = — 8 K 3 7 j(f), which leads 
to 

7 ,(t)= 7 #)e- 8 ^ = 7 ,e“ 8K?t - (4-20) 



(4.18) 

(4.19) 
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Then, from (2.9) and (4.15) we obtain the time evolution of the norming constants 
c\j(t) and c r j(t) as 

C\j(t) = c\j (0) e 4K ?‘ = eye 4 *? 4 , c Tj (t) = c Tj ( 0) e' 



” 4k j 1 = c r ,e 



-4 kU 



5. AKNS method to derive the KdV 

In the previous section we have seen that the KdV arises as a compatibility condi- 
tion in the Lax method. There are other methods to derive nonlinear PDEs that 
can be solved by the inverse scattering transform, i.e., by solving the inverse prob- 
lem with the time-evolved scattering data for a corresponding linear differential 
equation. One of these methods was developed by Ablowitz, Kaup, Newel, and 
Segur, and it was first applied to the Sine-Gordon equation [23]. Here we outline 
the basic idea behind the method of Ablowitz, Kaup, Newel, and Segur (AKNS 
method, for short) and use it to derive the KdV equation. 

Given a linear operator X associated with the system v x = Xv, we are 
interested in finding another operator T (the operators X and T are said to form 
an AKNS pair) such that: 

(i) The spectral parameter A does not change in time. 

(ii) The quantity v t — Tv must remain a solution to v x = Xv. 

(iii) The quantity X t -T x -{-XT — TX must be a (matrix) multiplication operator. 

Note that in general X contains the spectral parameter A, and hence T also de- 
pends on A as well. Usually, X and T are matrix- valued with entries depending on 
x,t, and A. As the operator A in the Lax method determines the time evolution of 
solutions to Cxjj = X'tp, in the AKNS method the operator T determines the time 
evolution of solutions to v x = Xv according to condition (ii) listed above. 

As the following argument shows, for compatibility, we are forced to have 

X t -T x + XT-TX = 0, (5.1) 

which leads to an integrable PDE and is in general nonlinear. From condition (ii) 
above we see that 

(v t - Tv) x = X(v t - Tv), 

or equivalently, 

vtx ~ T x v — Tv x = Xv t — XTv 

= (Xv) t - X t v - XTv 

V J (5.2) 

= (v x ) t - X t v - XTv V ; 

= v xt - X t v - XTv. 

We expect v to be smooth enough so that v tx = v xt . Let us replace Tv x by T Xv 
on the left-hand side in (5.2), from which we get (X t — T x + XT — T X)v = 0, 
which in turn, as a result of condition (iii) listed above, gives us the compatibility 
condition (5.1). 




12 



T. Aktosun 



Let us write the Schrodinger equation (1.2), by replacing the spectral param- 
eter k 2 by A, in the form of the first-order linear system v x = Xv by choosing 

|"0 u(x,t) — a] 



'Ipx 



x 



i 



0 



We will construct T so that T and X will form an AKNS pair. Let us try 

\a 0] 



T = 



£ V 



where the entries a , fi, £, and r) may depend on x, t, and A. The compatibility 
condition (5.1) leads to 



0 u t 

0 0 

or equivalently 



&x fix 

£ X Vx 



+ 



0 u — A 

1 0 



a fi 
€ V 



a fi 

£ V 



0 u — A 

1 0 



0 0 
0 0 



-a x - 0 + £(u - A) u t — 0 X + rj(u - A) - a(u - A) 




1 

O 

O 


+ a - T] -T) x + 0 - £(u - A) 




O ' 
O 



(5.3) 

From the (1, 1), (2, 1), and (2, 2) entries of the above matrix equation, we obtain 
fi = -a x + (u - A)£, rj = a-£ x , rjx = (5.4) 

Then the (1, 2) entry in (5.3) is given by 



Ut T 2 £>xxx 



u x £ - 2 6b ( u - A) = 0. (5.5) 

Letting £ = AC + \i in (5.5), where A is the spectral parameter, we obtain 

^CcA T ^~Cxxx ^CxU T 2 fi x u x (^ A T ^Ut T ~l^ xxx 2 n x u u x ia^ 0. 

Equating the coefficients of each power of A to zero, we get 
13 1 

C — c l, jU — 2 ClU + Ut - 2° lUUx “ + = 0, (5.6) 

where c\ and C2 are arbitrary constants. Using c\ = 4 and C2 = 0, from (5.6) we 
obtain the KdV given in (1.1). Moreover, with the help of (5.4) we get 

a = u x + c 3 , fi = -4A 2 + 2Au + 2ii 2 - u xa; , £ = 4A + 2u, rj = c 3 - 

where C3 is an arbitrary constant. Letting C3 = 0, we obtain 



T = 



— 4A 2 + 2A u + 2 u 2 



4A -(- 2ri 



u x 



(5.7) 



It is possible to obtain the time evolution of the Jost solutions and of the 
scattering data via the AKNS method. Note that condition (ii) of the AKNS 
method is satisfied if we have 



v t - Tv = p(t, A) v, 



(5.8) 
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for some (scalar) coefficient p(£, A). For example, if we choose 

_ f( (k, x; t) 

by letting x — > +oo in (5.8) and by using (2.2), we obtain the time evolution of 
f\(k,x\t) given in (4.13). By using (2.4) with T(x) replaced by T(k\t) and with 
L(k) by L{k\t ), and by letting x — > — oo in (5.8) and using (2.2), we get the time 
evolutions given in (4.15). In a similar way, with the help of (2.3), (4.9), and (5.7), 
by choosing 

I" fl(k,x-t) 
v = , 

fr(k,x;t)_ 

and by letting x — > ±oo in (5.8), we obtain (4.16) and (4.17). 

6. Solution to the Cauchy problem for the KdV 

The Cauchy problem (initial- value problem) for the KdV consists of finding u(: r, t) 
when u(x,0) is known. As shown by Gardner, Greene, Kruskal, and Miura [11], 
this problem can be solved by using the inverse scattering transform. Let V(t) 
denote any of the equivalents of the scattering data for the Schrodinger equation 
with the time-evolved potential u(x,t). In other words, we have 

m ■= {R(k;t), L(k;t),T(k]t), { 7j(0}> ( c y W}> ( c r j W)}> (6-1) 

or equivalently, 

V(t) = {R(k) e 8ik3t , L(k) e~ Sik3 \T(k), { Kj }, { 7j e^}, {cy e 4 ^}, { c rj e" 4 ^ 4 }}- 

Note that T>(0) corresponds to the initial scattering data associated with the po- 
tential V(x), where iz(x, 0) = V(x). 

Below we outline the solution to the Cauchy problem for the KdV. 

(i) Given u(x, 0), determine the corresponding scattering data Z>(0). This is done 
by solving the direct scattering problem V{x) i— ► T>( 0). The solution to this 
problem is essentially equivalent to solving (2.1) and obtaining the Jost so- 
lutions f\(k,x) and f r (k,x), from which 7^(0) can be constructed. 

(ii) Evolve in time the scattering data as V( 0) V(t) in accordance with (6.1). 

Note that the time evolution of the scattering data is really simple. On the 
other hand, the time evolution of the potential u(x, 0) i— ► u(x, t) will be much 
more complicated. Similarly, we expect that the time evolution of the Jost 
solutions f\(k,x) h-* f\(k,x\t) and f T (k,x) f r (k,x;t), governed by the 
PDEs (4.13) and (4.16), respectively, will be complicated. 

(iii) Having obtained the time-evolved scattering data V(t), solve the correspond- 
ing inverse scattering problem V(t) h-> tx(x, t) for (1.2). This problem is known 
to be uniquely solvable [13] when the initial potential V belongs to the Fad- 
deev class. 
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As indicated in Section 3, in step (iii) above one can obtain the solution to the 
inverse scattering problem by solving, for example, the time-evolved Riemann- 
Hilbert problem [cf. (3.1)] 

m\(—k, x ; t) = T(k) m r (fc, x\ t) — R(k) e 2lkx + 8lk 1 m\(k , x\ t), k £ R, (6.2) 
and recover u(x, t ) by using [cf. (3.4)] 



„ fet) = m; 'p ; ») +2 a fc€C+ , 



(6.3) 



m\{k,x\t) mi(/c,x;t)’ 

where the right-hand side can be evaluated at any k value, including k = 0 and 
k = dhoo. Equivalently, one can solve the time-evolved Riemann-Hilbert problem 
[cf. (3.2)] 

m T (-k, x; t) = T(k) rm(k, x\ t) - L(k) e - 2ikx - 8ikH m r (k, x-,t), k € R, 
and use [cf. (3.4)] 

= <(fc,g;t) _ m[(fc,x;f) £+ 

’ m r (k,x\t) m,;(k,x;t)’ ’ 

where the right-hand side can be evaluated at any k value in C+. 

Alternatively, one can solve the time-evolved left Marchenko equation [cf. 

(3.5)] 



with 



poo 

(x,y;t) + £l\(2x + y,t)+ dz£l](2x + y + z;t) B\(x,z\t) = 0, y > 0, (6.4) 

Jo 



1 poo iy 

dk R(k) e uk3t+ik y + Y, cl e' 

^ J — OO j = i 

and recover u(x,t ) by using [cf. (3.6)] 

dB\(x,0+]t) 



8 Kjt-Kjy 



u(x, t) = — 2 ■ 



<9x 



(6.5) 



Equivalently, one can solve the time-evolved right Marchenko equation [cf. (3.8)] 

poo 

B r (x,y;t) +tt r (-2x + y;t) + / cfo fi r (-2x + y + z; t) B r (x, z; t) = 0, y > 0, 

Jo 

( 6 . 6 ) 

with 

AT 

1 z* 00 

fi r (a;i):=— / dk L(k) e~ SikH - 

^ J —OO 

and obtain u(x,t) by using [cf. (3.9)] 

( dB r (x,0 + ;t) 






2 -8 Kjt-Kjy 

C r j e 



7 — 1 



u(x, t) = 2 



dx 



(6.7) 
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7. Soliton solutions to the KdV 



Consider the Cauchy problem for the KdV corresponding to the initial scattering 
data with zero reflection coefficients, N bound states at k = Ik j, and dependency 
constant 7 j, where we have the ordering 0 < K\ < • • • < kn- When R = 0, from 
(2.11) we see that the transmission coefficient is given by 



m = n 



3 = 1 



k + iKj 
k — iKj 



(7.1) 



In this case, the inverse scattering problem V(t) i-» u(x,t) can be solved alge- 
braically in a closed form, and the resulting solution u(x, t ) to the KdV is known 
as the V-soliton solution. Using (7.1) we can write (6.2) as 



N N 

m\(—k,x\t) fj(& — ify) = m r (k,x\t) JJ(& + iKj), k G R. (7.2) 

3 — 1 3 — 1 

Prom the analyticity properties of m\ and m r it follows that each side in (7.2) is 
entire in k with a polynomial growth of leading term k N as k — ► 00 in the complex 
plane C. Further, since k appears as ik in the Faddeev functions, both sides in 
(7.2) must be a polynomial of the form 

k N + ik N ~ 1 aN-i(x,t) + • • • + i N ao(x,t ), (7.3) 



where aj(x,t) are real valued and to be determined by using [cf. (2.10), 
(4.20)] 

e 2K jX -%K*t = j = l,...,N. 

m r (iKj,x ; t) 



(3.3), and 



(7.4) 



From (7.2) and (7.3) we get 



m\(k,x;t) 



k N — ik N 1 cln- i(x, t) + • • • + (— i) N ao(x, t) 

nf=i(fc + iK j) 



(7.5) 



m T (k,x;t) 



k N + ik N 1 cln-i ( x,t) + b i N ao(x,t) 



HL(k + i Kj ) 



Let us define 



Uj := e 2 ^ 1 " 8 ^*, 



j = 1,...,N, 



(7.6) 

(7.7) 



where we recall that the sign of ujj is the same [cf. (2.10)] as that of (— l) N ~i . 
Using (7.5)-(7.7) in (7.4) we get a system of linear algebraic equations for the N 
unknowns dj(x,t ), namely 



~ K ? 1 aN-i(x,t )H b (~l) N a 0 (x, t) 

+ K^ _1 ajv_i(a:, t) H b a 0 (x,t) 
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which can be written as 



r „n- 1 



(<* + 1) "V - 1) 
+ 1 ) ^ 2~ 2 ( U 2 - 1 ) 



ci-c-ir 

_l)iv 



U>2 



L^N 



N 1 (^N + 1) ftjv 2 ( lu N — 1) 



wn - (-1) 



N 





CljY—l (x, t) 




dN-2(x,t) 




a 0 (x,t) 



(7.8) 



r <(l-o;i) 

*2 (! “ ^2) 

L*n(1 — u n)\ 

Let M( n ^ denote the coefficient matrix in (7.8) whose (j, n) entry is given as 
Mjf := fa - (-1)”] , j, n = 1, . . . , TV. (7.9) 

Prom (7.8), via Cramer’s rule, it is possible to extract ajv_j(x,t) explicitly as 

det p( N ~i) 



a N -j(x,t) = — 



j = 1, . . . , TV, 



det MW ’ 

where the (TV + 1) x (TV + 1) matrix p( N ~ti is given by 

pW-i). = \ 0 Y(NJ) 1 

^ ■ |_QW mwJ ’ 

with <9< w ) being the TV x 1 matrix whose jth row contains the entry (1 —u>j) for 
j = 1, . . . , TV; and M (N> being the TV x TV matrix given in (7.9); and Y <N ' J> being 
the 1 x iV matrix whose jth column contains the entry 1 and all the remaining 



entries being zero. For example, we have 








0 


1 


0 


0 




(! - Wi) 


1 ( a, l + 1) 


Acf- 2 (^! - 1) ••• 


"i - (-1)" 


p(N- 1) . = 


K-2 (1 - w 2 ) 


K2~ l {u>2 + 1) 


^->2-1) ••• 


W2-(-l)" 




.Kjv(1 ~U N ) 


+ 1) 


k1 n~ 2 {un — 1) • • • 


UJn — ( — 1) W . 


We note that u(x,t) 


can be constructed from a at- i(#, t) alone 


: or from ao(x, 



alone. For example, using (6.3) and (7.5) in the limit k — > d=oo, we see that 

da N -i(x,t) 



u(x,t) = 2 ■ 



dx 



(7.10) 



or using (6.3) and (7.5) at k = 0 we see that 

1 d 2 ao(x,t) 



u(x,t) = 



ao(x,t) dx 2 



(7.11) 
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In fact, using (7.5) in (6.3) we see that one can recover u(x , t) by using any one of 
the following (N + 1) equations: 



u(x , t) = 



dN-j(x,t) L dx 2 



d 2 a N -j( x,t) da N -j-i^c,ty 



dx 



j = 0,1,..., AT, 



where we have defined a_i(x,t) := 0 and ajv(x,£) := 1. 

Alternatively, we can obtain the AT-soliton solution u(x, t) via (6.5) by solving, 
in a closed form, the left Marchenko equation (6.4), which has a degenerate kernel 
thanks to the fact that R = 0. Similarly, u(x, t) can be obtained via (6.7) by solving, 
in a closed form, the right Marchenko equation (6.6), which has a degenerate kernel. 

It is also possible to obtain the solutions to the Marchenko equations (6.4) 
and (6.6) in an algebraic manner, without really solving the integral equations 
themselves. Below we illustrate this for the recovery of B\(x, y\ t) and obtain u(x, t) 
via (6.5). Using R = 0 we can write (6.2) as 



m\(— fc, x; t) — 1 = [T(fc) — 1] m r (fe, x; t) + [m r (fc, x; t) — 1], k G R. 

Taking the Fourier transform of both sides with ^ d/c e lky and exploiting the 
analyticity properties of m\ and m r , we obtain [cf. (3.7), (3.10)] 

1 f°° 

B\(x,y;t) = — I dk[T{k)-l}m r (k,x-,t)e iky , y > 0. (7.12) 

^ J — on 



The right-hand side of (7.12) can be evaluated as a residue integral along the 
semicircular arc that is the boundary of C + . Using (7.1) in (7.12), we get 

N 

B\(x, y\ t) = i ^ [Res (T, iKj)] m T (iKj, x; t) e~ KjV , y > 0. (7-13) 

3 = i 

Using (2.9) and (7.4) in (7.13), we obtain 



N 

Bi(x,y;t) = -^^mi(i^,x;^)e" Kj ' y , (7.14) 

3 = 1 



where we have defined 

P . — r 2 -2 kj-x+ 8 K*t 

c 3 C 1 

Using the time-evolved version of (3.7), from (7.14) we get 



N 

m\(k,x-,t) = 1 — ^ - — pd — -m\(iKj,x]t) e~ Kjy . 
3=1 3 



(7.15) 
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Putting k = iK n in (7.15) forn = 1 , . . . , IV, it is possible to recover the m\(iK,j,x] t) 
by solving the linear algebraic system 



1 + 



gj 

«l + K 1 
gl 



^2 + «1 



1 + 



^2 
K 1 + 

^2 



^2 + « 2 



gj 

+ «1 



g2 

tt/V + «2 



£N 


1 




K 1 + I^N 




m\(iKi,x; t) 




T 


£N 




m\(iK2,x; t ) 




1 


K>2 + KN 






= 




1 + £N 




Tn\{in at , x\ t ) 




_ 1 _ 


&N + &N . 


1 





(7.16) 



Let T denote the coefficient matrix in (7.16), namely, let its (n, j) entry be given 
by 



where 5 n j denotes the Kronecker delta. Then, using (7.14) and (7.16), we obtain 
B\(x, 0 + ; t) as the ratio of two determinants as 



Bi(x,0 + ;t) = 



det Z 
detT’ 



where Z is the matrix defined as 

ro si 



Z:= 



1 1 + 



gi 



^2 

^2 



Kl + K 1 
gl 



«2 + « 1 



1 + 



«1 + ^2 
g2 



1 



gl 



K2 + ^2 



g2 



gjV 

gjV 



«i + ^at 
gjV 
+ tt/V 



1 + 



£n 



K>N + K 1 ACiV + ^2 

It can be shown that det Z = d x (det T) and hence with the help of (6.5), we get 



u(x, t) 



d_ 

dx 



1 <9(detT) 
det T dx 



(7.17) 



Some Mathematica notebooks based on ( 7 . 10 ), (7.11), and (7.17) for the 
evaluation of TV-soliton solutions to the KdV and their animations are available at 
the author’s web page [24]. 

We can use (7.10), (7.11), or (7.17) to analyze properties of solitons of the 
KdV. For example, when N = 1 , from (7.10) or (7.11) we get the single soliton 
solution to ( 1 . 1 ) as 

u(x, t) = — 2 k\ sech 2 (k\x — 4 n\t + y/ln 71 



(7.18) 
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It is seen that the amplitude of this wave is , it moves in the positive x direction 
with speed Ak\ , and the dependency constant 71 plays a role in the initial location 
of the soliton. The width of the soliton is inversely proportional to Kq, which can 
be seen, e.g., by using the fact that f^° dx yj —u(x, t) is equal to y/2ir. 

By exploiting the properties of one-soliton solutions to the KdV, one can 
show that as t — ► +00, the TV-soliton solution to the KdV resembles a train of 
N separate solitons each traveling with speed Ak!j. In this case, the KdV can be 
considered for all t e R and it can be shown that each soliton emerges from the 
nonlinear interaction by experiencing only a change in the phase. For details, the 
reader is referred to [3,4,6,25]. 

8. Conclusion 

In this section we will comment on three aspects of the KdV; namely, the time 
evolution stated in condition (ii) of the Lax method outlined in Section 4, the 
conserved quantities, and the Backlund transformation. 

Some references carelessly state the time evolution associated with the Lax 
method. For example, in (1.2.10) of [3, page 6] it is stated that the evolution of the 
solutions to (1.2) is given by d t ip = Aip with A as in (4.7), instead of the correct 
statement (ii) of Section 4. The incorrectness of d t ip = Aip can be demonstrated 
explicitly by an elementary example. Consider the 1-soliton solution to the KdV 
given in (7.18). Let us choose 71 = m = 1. A solution to (1.2) is obtained as 
f\(i,x]t). Let us call that solution ip. We have 

ip(x,t) = ^ e _4t sech(x — At), u(x,t) = — 2sech 2 (x — At). 

It can directly be verified that ip satisfies (1.2) with k — i, and (d t — A)ip = — Aip , 
as indicated by the correct time evolvement (4.13) with k = i; hence dtip A Aip. 

In (1.7.6) of [6, page 25] and (4) of [5, page 81] it is carelessly stated that 
the evolution of the solutions to (1.2) is given by dtip — Aip = cip with A as in 
(4.7) and c is an arbitrary constant or a function of t. The incorrectness of this 
can be demonstrated, for example, by choosing ip as f\(k,x;t) + f T (k,x;t), i.e., 
the sum of the Jost solutions to (1.2). From (4.13) and (4.16) it follows that, 
with this choice of ip , the equation dtip — Aip = cip would hold if and only if 
(Aik 3 + c)f r (k,x]t) — (Aik 3 - c)f\(k,x\t) = 0 for x G R and t > 0, which is 
impossible due to the linear independence of f\(k,x;t) and f r (k,x;t) on x G R. 
Note that (4.4) has two linearly independent (generalized) eigenfunctions for each 
A > 0 and hence dtip — Aip in general is not expected to be a constant multiple of ip. 

Let us note that a potential in the Faddeev class need not even be continuous. 
On the other hand, from (1.1) we see that classical solutions to the KdV are 
thrice differentiable with respect to x. Informally speaking, the discontinuities 
that may be present in the initial value u(x,0) disappear and u(x,t) becomes 
smoother for t > 0. On the other hand, even though u(x,t) changes as t increases, 
certain integrals involving u(x,t) with respect to x remain unchanged in time. 
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Such quantities are known as conserved quantities for the KdV. They can either 
be obtained directly from (1.1) or from the expansion of T(k) in powers of 1/ A; as 
k — > ±oo by using the fact that T(fc; t) does not change in time [cf. (4.15)]. When 
u(x,t) is smooth, with the help of (2.6) we obtain 



TIM -n a c? c? 

+ 2i k 8fe 2 4S.(" 

where we have defined 



+ ^ 3 +0(1/fc4) ’ 



k — > oo in C+, 



/ oo 

dxu(x,ty , j — 1,2. 

-oo 



Thus, we have identified two of the infinite number of conserved quantities; i.e., 
Ci and C 2 are independent of time and are equal to their values at t = 0. The 
time independence of C\ and C 2 can also be obtained directly from (1.1). We can 
write (1.1) as u t = (3 u 2 — u xx ) x , and hence 



dCi 

dt 



d_ 

dt 




dxu(x,t) = 3 u{x,t) 2 



'U'Xx(%-> t) 



00 

= 0, 

x= — oo 



where we have used u(x,t ) = o( 1) and u xx (x,t) = o( 1) as x — > ± 00 . Similarly, 
after multiplying (1.1) with u(x,t), we can write the resulting equation as 

(u 2 ) t = (4 u 3 - 2 uu xx + u 2 x ) x . (8.1) 



Integrating both sides of (8.1) on x G R, we get 

dC 2 
dt 



d f°° 

= — / dxu(x,t) 2 = 4u(x,t) 3 — 2u(x,t) u xx (x,t) +u x {pc,t) 2 
dt J —00 



= 0, 



verifying the time independence of C 2 . 

Can we characterize the set of nonlinear PDEs solvable by an inverse scat- 
tering transform? In other words, can we find a set of necessary and sufficient 
conditions that guarantee an initial-value problem for a nonlinear PDE to be solv- 
able via an inverse scattering transform related to a linear problem? There does 
not yet seem to be a satisfactory solution to this characterization problem. On the 
other hand, nonlinear PDEs solvable by an inverse scattering transform seem to 
have some common characteristic features such as the Lax pair, the AKNS pair, 
soliton solutions, an infinite number of conserved quantities, a Hamiltonian for- 
malism, the Painleve property, and the Backlund transformation. Here, we only 
briefly explain the last feature and refer the reader to [6] for details and other 
features. 

A Backlund transformation is a means to produce another integrable non- 
linear PDE from a given one. The basic idea is as follows. Assume v satisfies the 
integrable nonlinear PDE M{y) = 0, and u satisfies another integrable nonlinear 
PDE, say Q(u ) = 0. A relationship V(u,v) = 0, which is called a Backlund trans- 
formation, involving v, u , and their derivatives allows us to obtain Q{u ) = 0 from 
M.{y) = 0. A Backlund transformation can also be used on the same nonlinear 
PDE to produce another solution from a given solution. 
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As an example, assume that v satisfies the modified KdV given by 

v t - 6v 2 v x + v xxx = 0, x G R, t > 0. (8.2) 

Then, choosing 

u = v x + v 2 , (8.3) 

one can show that 

tit 6uu x T u xxx — ( 0 X T 2 V s ) (vt 6u v x T v xxx f x G R, t 0. 

Thus, (8.2) and (8.3) imply (1.1). The Backhand transformation given in (8.3) is 
known as Miura’s transformation [26] . For a Backlund transformation applied on 
the KdV to produce other solutions from a given solution, we refer the reader to 
[6,27]. 

Another interesting question is the determination of the linear problem asso- 
ciated with the inverse scattering transform. In other words, given a nonlinear PDE 
that is known to be solvable by an inverse scattering transform, can we determine 
the corresponding linear problem? There does not yet seem to be a completely 
satisfactory answer to this question. We mention that Wahlquist and Estabrook 
[28] developed the so-called prolongation method to derive the linear scattering 
problem associated with the KdV and refer the reader to [6] for details. 
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The Schur Algorithm for Generalized Schur 
Functions IV: Unitary Realizations 

D. Alpay, T.Ya. Azizov, A. Dijksma, H. Langer, and G. Wanjala 

Abstract. The generalized Schur transform as defined in [13] (see also [2]- 
[6]) is applied to the class A 0 of all complex- valued functions, which are 
holomorphic at z — 0. Each such function has a coisometric and a unitary 
realization in some Krein space. We study the effect of this generalized Schur 
transform to the unitary realization; in [2], [3] we studied similar questions for 
the coisometric realizations. The main difference with [2], [3] is that a certain 
one-sidedness is replaced by a two-sidedness, comparable to the difference 
between the unilateral shift on one-sided sequences and the shift on two-sided 
sequences. We follow a direct approach in line with [2, 3, 6]. 



1. Introduction 



1. Realizations. We denote by A 0 the class of functions which are holomorphic in 
a neighborhood of 0 in the open unit disk ID. Each function s(z) £ A 0 admits a 
unitary realization , that is, it can be written in the form 

s(z) = 7 + z((Ijc - zT)~ 1 u, v) K . (1.1) 

Here K is a Krein space with indefinite inner product ( • , • )*;, Ijc is the identity 
operator and T is a bounded operator on /C, u and v are elements in /C, and 7 is 
a complex number: 7 = s(0), such that the matrix operator 



( T u 
\( * > v )/c 7, 




X 

C 



X' 

C 



/c®c, 



(1.2) 



is unitary: U* = U~ 1 . We recall the terminology: The matrix U is identified with 
the colligation {/C, C, T, u , v, 7}, X is called the state space , T is the main operator 
and the function on the right-hand side of (1.1) is called the characteristic function 
of the colligation. Evidently, the characteristic function of a unitary colligation 
belongs to the class A 0 . 
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That the colligation is unitary can be expressed by the relations 



TT* + (-,u) k u = I K , (1.3) 

(v,v)k + \j \ 2 = 1, (1.4) 

Tv + i*u = 0, (1.5) 

T*T+(-,v) k v = Ik, (1.6) 

( u,u) K + \-y\ 2 = 1, (1.7) 

T*u + jv = 0. (1.8) 



The unitary colligation (1.2) can be chosen minimal which means that 

sp n{T n u, T™v | n, m = 0, 1, 2, . . .} = /C. 

If this holds then (1.1) is called a minimal realization of s(z). A minimal unitary 
realization of a function s(z) £ A 0 is essentially unique: If besides (1.1) also 

s(z) = 7 + z ((/~ - zfy^v)^. 

is a minimal unitary realization of s(z), then there is a weak isomorphism W from 
/C to /C (by definition, an operator mapping a dense linear subset of /C isometrically 
onto a dense linear subset of JC) such that 

WT Idom w = TW , Wu = 5, Wv = V. 

In this case we have that if /C = /C+ 0 /C_ and /C = /C+ 0 /C_ are any fundamental 
decompositions of K and /C respectively, then dim 1C± = dim K ± , in the sense that 
either both numbers are infinite or both are finite and equal. If at least one of 
these numbers is finite then W isJ>ounded and can be extended by continuity to 
a unitary mapping from /C onto /C and WT = TW. These facts follow from, for 
example, [10, Section 5.3] and [14, Theorem 2.2]. We assume the reader is familiar 
with Krein and Pontryagin spaces and the theory of operators on these spaces such 
as in [10], [17], and [15]. 

2. Generalized Schur functions. Of the many equivalent definitions of a generalized 
Schur function from the class S° (see [18] and also [7]) we select the one directly 
related to the preceeding subsection: s(z) belongs to the class S° if s(z) £ A 0 and 
the state space in the minimal unitary realization of s(z) is a Pontryagin space, 
that is, a Krein space /C for which, in the notation above, dim /C_ < oo. The 
subclass of functions s(z) £ S° for which dim /C_ = k < oo, will be denoted by 
S^. If k = 0, the state space is a Hilbert space and s(z) belongs to Sq if and only 
if it can be extended to a Schur function, that is, a holomorphic function defined 
on all of B which is bounded by 1; we write So for S§. Finally, we mention that 
s(z) £ if and only if it admits the representation 

s (*) = fl 

j = 1 3 

with points Zj € B, z 3 ^ 0, and a function so(z) £ Sq, so(zj) 0, j = 1, 2, . . . , k. 
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3. The generalized Schur transform. In the sequel we write the Taylor expansion 
of s(z) G A 0 at z = 0 as 

s(z) = a 0 -ha 1 z-j-a 2 z 2 -i h cr n z n -i , a n = s (n) (0)/n!, n = 0,1, (1.9) 



The generalized Schur transform s(z) of a function s(z) is defined via a fractional 
linear transformation which depends on whether |<Jo| < 1 , |oo| > 1 , or |ctq | = 1: 






1 s(z) - <T 0 
Z 1 — CTqS(z) 

t : .k l ~ °*o s ( z ) 

s(z)-a o 

q (Q(z) - Z k )s{z) - Q-qQ(z) 
k #0 Q(z)s(z) - {Q(z) + z k ) 



if \a 0 \ < 1, 
if | <7 0 1 > 1, 
if | <7o I = 1 



( 1 . 10 ) 



and this equality holds for z in a neighborhood of 0 in D. The Schur transform is 
defined for those functions s(z) G A 0 which satisfy 

(i) s(z) ^ cr 0 with |cr 0 | > 1, and 

(ii) °oQ(z)s(z) - ( Q(z ) + z k )=£ 0 if |<r 0 | = 1. 

In this definition k is the smallest integer > 1 for which o k ^ 0. The condition 
(i) implies that such a k exists. The polynomial Q(z) of degree 2k and the integer 
q > 0 in case fyo| = 1 are defined as follows. Introduce the complex numbers c n , 
n = 0, 1, . . ., by the relation 



(s(z) — <Jo)(c 0 + C\Z + * * • + C n Z n + •••)= CFQZ k , 



then Q(z) is defined by 



Q(z) = P(z) - z 2k P(l/z*)*, P(z) := Co + c\z H V c k -\z k 1 . 

It is easy to see that (Q(z) — z k )s (z) — ctqQ(z) has a zero at z = 0 of order 2k and 
that ctqQ(z)s(z) — ( Q(z ) + z k ) has a zero at z = 0 of order at least 2k. On account 
of condition (ii) the quotient of these functions has a pole at z = 0 of finite order. 
In the definition of the generalized Schur transform (1.10) this order is denoted 
by q. 

It can be shown (see [16], [11], and also [2]) that a generalized Schur function 
which is not identically equal to a unimodular constant satisfies (i) and (ii) and 
that the generalized Schur transform of such a function is a generalized Schur 
function. 

The Schur transform goes back to [19] and [20] in connection with the study 
of interpolation and moment problems. The generalized Schur transform was first 
developed by the authors of [12] and [16] in order to generate certain Pisot num- 
bers. Later the transform was perfected to the present definition; for details see 
[11]. In [13] it was studied in a system theoretic context. The authors of [8] and 
[9] also investigated the Schur transform in the nonpositive case. They considered 
matrix- valued functions (in this connection see also [1]) and used reproducing 
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kernel methods. They “avoided” the case where |<7o| = 1 by applying the Schur 
transform in points of D different from the origin. 

4. Contents of the paper. In this paper we investigate how to obtain the minimal 
unitary colligation of the Schur transform s(z) of s(z) from the minimal unitary 
colligation of s(z). In Section 2 we formulate the theorems: Theorem 2.1 for the 
case |cr 0 | < 1, Theorem 2.3 for the case |cr 0 1 > 1, and Theorems 2.4 and 2.5 for 
the case |<Jo| = 1, the first one when q = 0 and the second one when q ^ 0. The 
final Theorem 2.6 concerns the minimal unitary realization of the function which 
is obtained by applying the generalized Schur transform finitely many times to 
5 ( 2 ). In Section 3 we provide proofs, for Theorems 2.1 and 2.3 rather briefly, for 
the Theorems 2.4 and 2.5 in detail. 

The results are similar to the results in [2] and [3] where we studied the 
effect of the Schur transform on the minimal coisometric realizations : s(z) has 
the representation (1.1) but now the colligation U in (1.2) is (i) coisometric , that 
is, UU* = I jc or, equivalently, the relations (1.3), (1.4), and (1.5) hold, and (ii) 
minimal , that is, 

span {T* m v \ m = 0, 1, 2, . . .} = 1C. 

In [2], [3] a model for the main operator in the minimal coisometric colligation 
was obtained by considering the sequence 

v , T*v, T* 2 v,...,T* m v,... 

as a basis for /C. The elements ^2 &nT* n v in the span of this basis were identi- 
fied with the one-sided sequences (a n )^T 0 with finitely many nonzero terms; for 
example, v was identified with (1, 0, 0, . . .). Then T* is just the right shift on the se- 
quence space. The inner product on /C was expressed in terms of the inner product 
on the space £ 2 : 

{^a n T* n v,Y,!3nT m v)^ = (0(a„)S°,(/? n )8% , 

where the Gram matrix is given by Q = X — T T* with indentity matrix X and 





o- 0 


0 


0 • 


•• \ 




01 


00 


0 • 






02 


01 


a 0 ■ 












■ / 



In this note as a basis for the state space of the minimal unitary colligation U we 
choose the two sided sequence 

...,T* k v,...,T* 2 v,T*v,v | u,Tu,T 2 u,...,T e u,.... 

The elements in its span can be identified with two-sided sequences with finitely 
many nonzero terms. For example, 

« = (■•■ ,0,0 1 1,0- - « = (••• ,0, 1 1 0,0, • •). 
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The model for the main operator T is the matrix 



1 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 ••• 


~ a 2 




~ a 0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


1 



where <to = s(0) and Oj+\ = (T j u, v)jc is the coefficient of in the Taylor series 
expansion of s(z), j = 0, 1, 2, . . .. By the relations (1.3)-(1.8) (see Lemma 3.1 in 
Section 3), the Gram matrix corresponding to this two-sided sequence is given by 




’ I -AA* 


B 


c 


1 - V*V 



with the infinite Hankel matrices 



/, : : :\ 





<?0 


01 


<?2 






0 


00 


o \ 




v- 


0 


0 


<70/ 




/... 


_* 

a 3 


a 2 




\ 




(T4* 


a 3 


a 2 




*5 


a 4 


a 3 



V ■ ■ ■ ) 



( : : : \ 



03 


04 


05 * 




02 


03 


04 • 




\ 


02 


03 ‘ 


• / 




(<70 


01 


02 


...\ 




0 


00 


01 






0 


0 


00 















Thus, unlike in the coisometric model, in the unitary model the roles of the ele- 
ments u and v are similar. This is reflected in the theorems that follow. 



2. The theorems 



In the following theorems we consider a function s(z) £ A 0 with Taylor expansion 
(1.9): 

s(z) = 00 + <j\z + a 2 z 2 -| b a n z n H , a n = s (n) ( 0)/n!, n = 0, 1, ... , 



U = 



(2.1) 



and assume 

T u\ . (K\ (K> 

K {-’ v )>c 7/ ' ^ vc, 

is the minimal unitary colligation, whose characteristic function coincides with 
s(z), that is, in a neighborhood of z = 0 in D we have 

s(z) = 7 + z((I , c - zT)~ l u, v)k, 7 = s(0) = a 0 . 



Here the state space K of U is a Krein space. 
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1. Case | <jo | < 1. In this case, let & E { 1,2,.. .} and <ti ■ 
Unlike in the next case we do not require that Gk ^ 
transform k times to s(z) then we obtain the function 

1 s(z ) - Gq 



s(z) 



Z k 1 - CTqS(z) ' 



= • • • = Gk-i = 0 if k > 2 . 
0. If we apply the Schur 



(2.2) 



We express the colligation of this function in terms of that of s(z). 



Theorem 2.1. Assume s(z) E A 0 and \go\ < 1. Let k be an integer > 1 such that 
G\ — ... = g k—i = 0 if k > 2. Let U in (2.1) be the minimal unitary colligation 
whose characteristic function coincides with s(z). Then 

C = span {v, T*v, T^ k ~^v}, £ = span {u, Tu , . . . , T (A:_ 1 ) u}, 

are k-dimensional positive subspaces of JC and s(z) defined by ( 2 . 2 ) is the charac- 
teristic function of the minimal unitary colligations U and U : 



with 



T 

v 



U = 



T u\ (K\ (K 

7/AC7 VC 



PTP, 



V 1 - kol 



1 PT* k v, 7 



. — Pu , 



where P is the orthogonal projection in JC onto the subspace JC := JC © C, and 



U = 



u\ (JC 



( * > v )jr 7/ 



with 



QTQ, 



u — 



-Qv, 7 = 



V 1 - Kl 2 



QT k u, 



' i-koi 2 ’ 

where Q is the orthogonal projection in JC onto the subspace JC := JC © C. 



Remark 2.2. The linear map W from 

span {T n 2, T* m v | n, m = 0, 1, 2, . . .} C JC 

onto 

span {T n u, T^v \ n, m = 0, 1, 2 , . . .} C JC 
defined by WT n u — T n u , WT™v = T* m 7, n, m = 0, 1, 2, . . ., is a weak isomor- 
phism from JCtoJC such that 

wf Idom w= TW, Wu = 2, Wv = V. 

This remark also applies to Theorems 2.3 and 2.5 below. A proof is given at the 
end of Section 3. 
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The theorem implies that the state space /C loses k positive squares. Thus if 
s(z) belongs to and satisfies the conditions of the theorem, then s(z) defined 
by (2.2) also belongs to S°. 



2. Case |<7o| > 1. We assume s(z ) is not identically equal to ao- Then (and only 
then) there is a smallest integer k > 1 such that ak / 0. Thus k E {1,2,...} and 
(a) G\ ^ 0 if k = 1 and (b) <ti = . . . &k-i = 0 and <Jk 7 ^ 0 if k > 2. The Schur 
transform of s(z) is given by 



s(z) = z k 



1-^0 s(z) 
s(z) - a 0 ' 



We express the colligation of s(z) in terms of that of s(z 



(2.3) 



Theorem 2.3. Assume s(z) E A 0 , |cr 0 1 > 1 , and s(z) ^ <7o. Let k be the smallest 
integer > 1 such that ^ 0. Let U in (2.1) be the minimal unitary colligation 
whose characteristic function coincides with s(z). Then 

C = span {i?, T*v , . . . , T*^ k ~^v}, C = span { u , Tu , . . . , T^ k ~^u}, 

are k-dimensional negative subspaces of 1C and's(z) defined by (2.3) is the charac- 
teristic function of the minimal unitary colligations U and U : 



U = 



u\ 1C 



V 'Vc 



with 



T = PTP — 



VKI 2 ~ 1 



ffc 

~-PT* k v, 



Pu , u 



7 = 



VWo \ 2 - 1 

& k 

1 - Kl 2 



Pu , 



a k ak 

where P is the orthogonal projection in 1C onto the subspace 1C := ICQ C, and 

u\ (K\ (K> 



u = 



T 

( • i v )jc 1) \c 



with 



QTQ- 



{ - ,Qv) ^ QT k u, u 
<*k 



V = — Qv, 



7 = 



v' I<, o |2 ~ 1 0 T‘ u, 

1-kol 2 , 

&k 



where Q is the orthogonal projection in 1C onto the subspace 1C := ICQ C. 



The theorem implies that the state space 1C loses k negative squares. Thus 
if s(z) belongs to S° and satisfies the conditions of the theorem, then 1 < k < k, 
and s(z) defined by (2.3) belongs to S° K _ k . 
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3. Case \do\ = 1. As in the previous case we assume that s(z) is not identically 
equal to do and denote by k the smallest positive integer such that dk ^ 0. Let 
Co, ci, . . . , c n , . . . be the complex numbers satisfying 

(s(z) — cr 0 ){co + C\z + • • • + c n z n + •••) = &oz k , (2-4) 

and set 

Q(z) := P(z) - z 2k P(l/z*)*, P{z)~c 0 + c 1 z + --- + c k - 1 z k - 1 . (2.5) 



Note that Co 7^ 0, so that P( 0) 7^ 0 and Q(z) is a polynomial of degree 2k. We also 
assume that s(z) is such that <Tq Q(z)s(z) — (Q(z) 4 - z k ) is not identically equal to 
0. Then the Schur transform of s(z) is defined by 



~ M = <, (Q(z) - Z k )s(z) - a 0 Q(z) 
a oQ( z )s(z) - ( Q(z ) + z k y 



(2.6) 



where q is the order of the pole of the quotient on the right-hand side; hence q 
is a nonnegative integer. At z — 0 the numerator has a zero of order 2k and the 
denominator has a zero of order 2k + q. We write the Taylor expansion at z = 0 
of the latter as 



a* 0 Q(z)s(z)-(Q(z) + z k ) = t 2k+q z 2k+ i + O(z 2k+ * +1 ), t 2k+q ^ 0. (2.7) 

In the following two theorems we express the colligation of s(z) in terms of that 
of s(z). The first theorem concerns the case q — 0, which means that 



t2k ~ (C0°2 k + Cl&2k-1 + * ' ' + Ck-\dk+l) 7^ 0. (2-8) 

The second theorem deals with the case q > 0. 



Theorem 2.4. Assume s(z) E A 0 , |<7o| = 1, and s(z) ^ cjq. Let k be the smallest 
integer > 1 such that dk / 0 and assume (2.8) holds. Let U in (2.1) be the minimal 
unitary colligation whose characteristic function coincides with s(z). Then 

C 0 = span . . . , T*v, v, u, Tu , . . . T k ~ l u] 



is a (2k) -dimensional subspace with k positive and k negative squares, and s(z ) 
defined by ( 2 . 6 ) with q = 0 is the characteristic function of the minimal unitary 
colligation 



with 



U = 



T u\ f K\ (1C 



T 



PTP 



■ , PT* k v) 
&kt2k 



-PT k u, u 



—PT k u, 

t2k 



V 



- PT* k v , 



L 2k 



00 — 7 — ^ 
t2k 



where P is the orthogonal projection in 1C onto the subspace 1C := 1C Q Co- 
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Theorem 2.5. Assume s(z) G A 0 , |cr 0 1 = 1, and s(z) ^ <Jo- Let k be the smallest 
integer > 1 such that ^ 0 and assume that for some integer q > 0 

°0 Q( z ) s ( z ) — ( Q( z ) + Z k ) = t2k+qZ 2k+q + t2k+q+\Z 2k+q + l + * * * , ^2/c+q 7 ^ 0. 

Te£ U in (2.1) be the minimal unitary colligation whose characteristic function 
coincides with s(z). Then 

C q = span {r*( fc +«- 1 ) w , . . . , T*V, v, u, Tu,... T k ~ l u} 

and 

Z q = span , T*v, v, u, Tu,... T k+q ~ l u } 

are (2k + q)~ dimensional subspaces of 1C with k positive and k + q negative squares, 
and s(z) defined by (2.6) is the characteristic function of the minimal unitary 
colligations U and U: 



with 



U 



FTP - 



T u 
<•>”)£ 7 



h k+q 

— 3 : pT*(k+q) 

l 2 k+q 



/ PT < 1 c +*) v ) ? . 

— —PT k u, u 



1 



t2 k+q 
<?k 



- PT k u , 



^2 k+q 



where P is the orthogonal projection in 1C onto the subspace 1C := 1C 0 C q , and 



T u\ . (K\ ^ (K 
(•»£ 7/ V c / V c 



with 



_ /. PT* k v) - 

T = QTQ - — S 

v = 7 

C 2/e+g 

where Q is the orthogonal projection in 1C onto the 



= — - — QT k+q u, 

^2 k+q 

— (Jk 
^2k+q 

subspace 1C := 1C © C q . 



The two subspaces 

span {T< k ~ l \, ..., T*v, v}, span {u, Tu,... T fc-1 «} 

of Co in Theorem 2.4 and of C q and C q in Theorem 2.5 are neutral and skewly 
linked; see Lemma 3.4 below. 

The theorems imply that the state space /C of U loses k positive and k + q 
negative squares. Thus if s(z) belongs to S° and the case |cr 0 | = 1 applies then 
1 < k + q < k and the Schur transform s(z) defined by (2.6) belongs to . 
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4. The generalized Schur algorithm. Assume s(z) G A 0 satisfies the conditions 
(see Subsection 1.3): 

(i) s(z) ^ <7o with | oo | > 1, and 

(ii) <JqQ(z)s(z) - ( Q(z ) + z k ) ^ 0 if |<r 0 | = 1. 

Repeated application of the Schur transform to s(z) leads to a sequence of func- 
tions: 

s 0 (z) :=s(z), si(z) ~s 0 (z),...,Sj(z) :=Sj-i(z),..., 

This is called the generalized Schur algorithm. The algorithm stops with the func- 
tion Sj(z) when this function no longer satisfies the conditions (i) and (ii). If s(z) 
belongs to the class S° then this is the case when Sj(z) is identically equal to a 
constant of modulus 1. 

Theorem 2.6. Assume s(z) G A 0 satisfies the conditions (i) and (ii) and let U 
in (2.1) be the minimal unitary colligation whose characteristic function coincides 
with s(z). Let 3(z) be one of the functions appearing in the Schur algorithm of 
s(z). Then there exist integers t, r > 0 and complex numbers a ^ 0,6 / 0, and c, 
which may be 0, such that 

C = span . . . , T*v, v, u, Tu, . . . T r_1/ u} 

is a Pontryagin subspace of K, and s(z) is the characteristic function of the min- 
imal unitary colligation 




with _ 

T = PTP-c(,PT H v) K PT r u , u = aPT r u , 

v = bPT* l v , 7 = 3(0), 

where P is the orthogonal projection in /C onto the subspace /C := /C © C. 

The proof of Theorem 2.6 is by induction on the steps of the generalized Schur 
algorithm. According to Theorems 2.1, 2.3, 2.4, and 2.5, this theorem is valid if 
s(z) is the immediate successor of s(z). Assume that s(z) is as in the theorem and 
that the Schur transform can be applied to it. Then by invoking the theorems just 
mentioned one can show that the Schur transform of s(z) has a realization similar 
to that of s(z). We omit the details. 

3. Proofs of the theorems and the remark 

In this section we give the proofs of the theorems in the previous section. Since 
a unitary colligation is also coisometric, in our proofs we shall frequently refer to 
the theorems and the calculations given in [2] in which the minimal coisometric 
versions of these theorems are treated. We shall give complete proofs of Theorems 
2.4 and 2.5 as they have no analogs in the coisometric setting. 
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Consider s(z) E A 0 and assume that it has the minimal unitary realization 
(1.1) and the Taylor expansion (1.9) around z = 0: 

s(z) = 7 + z((Ijc ~ zT)~ l u , v)jc = (Jo + &iz ~\ h ojz 3 H , 

whence cr 0 = 7 and = (T j u,v)jc^ j = 0,1,2, The relations (1.3)-(1.8) 

imply the following lemma. 

Lemma 3.1. For i,j = 0,1,2,... it holds that (T' l u,T* j v})c = (Ji+j+i and 
{T^Tiujjc = (T*iv,T*'v)jc = 

= f (1 - M 2 - |(Ji| 2 - ... - |<r*| 2 ), i=j> 0, 

\ ~( a o < 7 i-j + otai-j+i + • • • + o’jai), i > j > 0 . 

For the last equality in the lemma see [2, Formula (4.7)]). From this lemma 
one can derive the formulas for the Gram matrices in Subsection 1.4. 



1 . Proofs of Theorems 2.1 and 2.3. In the proof of Theorem 2.1 we use the following 
corollary of Lemma 3.1. 

Lemma 3.2. If <ti = . . . = dk-i = 0 then for k > 1, the Gram matrices associated 
with the finite sequences 

v,T*v,...,T*( k -Vv and u,Tu, . . . ,T k ~ l u 

in 1C coincide and are equal to (1 — \ao\ 2 )Ik, where Ik is the k x k identity matrix, 
and {T* k v,T* k v) K = {T k u,T k u) K = 1 - |< 7 0 | 2 - M 2 . 



Proof of Theorem 2.1. The statements concerning the colligation U follow from 
those concerning the colligation [/, by applying the latter to the function s*(z) := 
s(z*)*: s*(z) is the characteristic function of the minimal unitary colligation U* 
and the Schur transform of s*(z) is given by (s)*(z). (Note that the roles of u and 
v are interchanged.) 

As to the first part of the theorem we first observe that by Lemma 3.2, C is a 
^-dimensional positive subspace. That U is coisometric and that its characteristic 
function coincides with s(z) follows from [2, Theorem 5.1] for the case k = 1 and 
[2, Corollary 5.2] for the case k > 1, proved by induction on k. The relations (1.3) 
and (1.5) for U are proved explicitly there. The remaining relations can be shown 
in a similar way. We show this for (1.7): Using the same relation but then for U 
we obtain 



(u,u)jc + |7| 2 = — j—^{ Pu ,u) K + |7| 2 



1 

1 

i^RF 



(u,u)tc - 



\(v,u)k\ 



1-KI2- 



{v,v) K 

N 2 



i-KI 2 




= l. 
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Thus U is a unitary colligation. It remains to show its minimality. For that we 
first note that there are constants a m j and /3 m j such that for m = 0, 1, 2 , . . ., 

m— 1 m— 1 

(. PT*) m v = T* m v + ^2 (XmjT* k v, ( TP) m u = T m u + ^ /3 mj T k u. 

3 = 0 3=0 

This can be seen by induction and the use of (1.3)-(1.8). Thus for x G 1C we have 

m— 1 

(T* m v,x)£. = {T*(PT*) m v, x)k. = {T< m+ ^v,x) K + £ a mj {T*U +l \x) K 

3=0 

and 

m— 1 

{T m ux)£ = ((PT) m u,x)ic = {T m u,x) K + P mj {T j u,x) K . 

3 = 0 

Induction arguments again show that for x G /C, 

(f* m u,x)£ =0, m = 0,1,2,... => (T* m v,x)^ =0, m = 1,2,..., 

(f m u,x)^ = 0, m = 0, 1, 2, . . . => (T m w,x)x:=0, m = 0,l,2,.... 

These implications in turn imply that the minimality of U follows from that of U. 

□ 



In the proof of Theorem 2.3 we follow the proof of its coisometric analog [2, 
Theorem 6.1]. We use the following lemma; see [2, Lemma 4.3] for its proof and 
references to earlier results. By the anti-space of a Krein space (/C, ( • , • )) we mean 
the Krein space (1C', ( • , • )') where K' — K as linear spaces and (x, y) f = — (x, y), 
x, y G /C'. Evidently, (/C')' = JC as Krein spaces. 



Lemma 3.3. Let f(z) be the characteristic function of a minimal unitary colligation 
U on 1C ® C of the form (1.2) and assume 7 = /( 0) ^ 0. Then the reciprocal 
function f r (z) := f(z)- 1 is the characteristic function of the minimal unitary 
colligation 



U r = 



T r 

(-,v r y 



u r 1 
7 r ) 



on /C' © C, where 1C' is the anti- space of the Krein space 1C, with entries 



T 
± r 

V r 



T - U r 



7 



V 

/*y * 



u 

5 

7 

1 



that is, fr(z) =7 r + Z ((Ijc — zT r ) 1 U r , V r ) f . 



7i 



7 
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Proof of Theorem 2.3. We rewrite the formula (2.3) in the form 

1 _1_ 

1 _ J0_ ]_ S{z) (Tp 

I(z) 



<7n 2* 



1 



(3.1) 



1 s ( z ) 



1 z k + 0(z k+1 ), 



and note that 1 1 / cro | < 1 and that 

1 1 

s(z) (Jo ^0 

in particular, this difference has a zero at z = 0 of order k. We apply Lemma 3.3 to 
s(z) and the corresponding colligation U and then we apply Theorem 2.1 to l/s(z). 
The formula (3.1) implies that ofi / (oo’s(z)) is the generalized Schur transform of 
1 /s(z)\ see (2.2). If (/C',T, u,v,^) is the minimal colligation for Oq / (jos(z)) , then 

a 0 CTo Jo 

is a minimal colligation for 1/J(z). Finally we apply Lemma 3.3 to 1 /s(z) and this 
colligation and obtain the realizations mentioned in the theorem. □ 

2. Proof of Theorem 2.4. We begin with the following consequence of Lemma 3.1. 
For k complex numbers zi, Z 2 , . . . Zk we define the k x k Toeplitz matrix 

\ 



A(z 1 ,z 2 ,...,z k ) 



( 

0 



0 

V o 



22 

Zi 

0 

0 



Zk - 1 
Zk- 2 



z 1 

0 



Zk 

Zk- 1 

Z2 

Zl 



J 



Lemma 3.4. The following three statements are equivalent: 

(1) | (Jo | = 1 and G\ — . . . = Ok - 1 = 0. 

(2) span{T*^ -1 ^, . . . ,T*v,v} is a k-dimensional neutral subspace of 1C. 

(3) span {u,Tt/, . . . ,T fc_1 u} is a k-dimensional neutral subspace of 1C. 

If ( l)-(3) hold then the spaces in (2) and (3) are skewly linked if and only if Ok ^ 0. 
In this case, the Gram matrix Q for the sequence 






, T*v, v | u,Tu , . . . 



,T k - x u 



is the 2k x 2k matrix 

g = 



o 

A* 



A' 

0 



A = A((jfc,<jfc+i, . . . ,02fc-i); 



in particular, the space 

span , T*v, v, u,Tu,..., T fe_1 u} 

is a (2k) -dimensional Pontryagin space with positive index k and negative index k. 




